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We derive quantum kinetic equations from a quantum field theory implementing a diagrammatic 
perturbative expansion improved by a resummation via the dynamical renormalization group. The 
method begins by obtaining the equation of motion of the distribution function in perturbation the- 
ory. The solution of this equation of motion reveals secular terms that grow in time, the dynamical 
renormalization group resums these secular terms in real time and leads directly to the quantum 
kinetic equation. This method allows to include consistently medium effects via resummations akin 
to hard thermal loops but away from equilibrium. A close relationship between this approach and 
the renormalization group in Euclidean field theory is established. In particular, coarse graining, 
stationary solutions, relaxation time approximation and relaxation rates have a natural parallel as 
irrelevant operators, fixed points, linearization and stability exponents in the Euclidean renormal- 
ization group, respectively. We used this method to study the relaxation in a cool gas of pions and 
sigma mesons in the 0(4) chiral linear sigma model. We obtain in relaxation time approximation the 
pion and sigma meson relaxation rates. We also find that in large momentum limit emission and ab- 
sorption of massless pions result in threshold infrared divergence in sigma meson relaxation rate and 
lead to a crossover behavior in relaxation. We then study the relaxation of charged quasiparticles 
in scalar electrodynamics (SQED). We begin with a gauge invariant description of the distribution 
function and implement the hard thermal loop resummation for longitudinal and transverse photons 
as well as for the scalars. While longitudinal, Debye screened photons lead to purely exponential 
relaxation, transverse photons, only dynamically screened by Landau damping lead to anomalous 
(non-exponential) relaxation, thus leading to a crossover between two different relaxational regimes. 
We emphasize that infrared divergent damping rates are indicative of non-exponential relaxation 
and the dynamical renormalization group reveals the correct relaxation directly in real time. Fur- 
thermore the relaxational time scales for charged quasiparticles are similar to those found in QCD 
in a self-consistent HTL resummation. Finally we also show that this method provides a natural 
framework to interpret and resolve the issue of pinch singularities out of equilibrium and establish 
a direct correspondence between pinch singularities and secular terms in time-dependent perturba- 
tion theory. We argue that this method is particularly well suited to study quantum kinetics and 
transport in gauge theories. 



I. INTRODUCTION 



The search for the quark-gluon plasma (QGP) at the Relativistic Heavy Ion Colhder (RHIC) and the forthcoming 
Large Hadron Colhder (LHC) has the potential of providing clear evidence for the formation of a deconfined plasma 
of quarks and gluons and hopefully to study the chiral phase transition. Perhaps this is the only opportunity to 
study phase transitions that are conjectured to occur in particle physics with earth-bound accelerators and an intense 
theoretical effort has developed parallel to the experimental program that seeks to understand the signatures of the 
QGP and the chiral phase transition An important part of the program is to assess whether the plasma, 

once formed, achieves a state of thermodynamic equilibrium and if so on what time scales. This is an important 
question since current estimates suggest that at the energies and luminosities to be achieved at RHIC, the spatial and 
temporal scales for the existence of the QGP are of the order of 20 fm The description of the space-time evolution 
in an ultrarelativistic heavy ion collision requires the understanding of phenomena on different time and spatial scales. 
Ideally, such a description should begin from the parton distribution functions of the colliding nuclei as the initial state 
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and evolve this state in time using QCD to obtain the kinetic and chemical equilibration of partons, the emergence of 
hydrodynamics and the hadronization and freeze-out stages |^ . An important part of the program to study the space- 
time evolution from first principles seeks to establish a consistent kinetic description of transport phenomena in a dense 
partonic environment. Such kinetic description has the potential of providing a detailed understanding of collective 
flow, observables (hadronic and electromagnetic) such as multiparticle distributions, charmonium suppression, freeze 
out of hadrons and other important experimental signatures that will lead to an unambiguous determination of whether 
a QGP has been formed and the observables of phase transitions. This premise justifies an important theoretical effort 
to obtain such a kinetic description from first principles. During the last few years there have been important advances 
in this program, from derivations of kinetic and transport equations from first principles in QCD and scalar 

field theories |]7|-[l0| to numerical codes that describe the space-time evolution in terms of partonic cascades [|| that 
include screening corrections in the scattering cross sections |ri|,p^ and more recently nonequilibrium dynamics has 
been studied via lattice simulations 13|-1^ 



The kinetic description to study hot and/or dense quantum field theory systems is also of fundamental importance 
in the understanding of the emergence of hydrodynamics in the long- wavelength limit of a quantum field theory 
and more recently a transport appro ach has been advocated as a description of the collective dynamics of soft degrees 
of freedom in hot QCD [p^|Jlqj20|,|2^,pl[| . The typical approach to derive transport equations begins by introducing a 
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Wigner transform of a particular nonequilibrium Green's functions at two different space-time points y-|5|,|20|,E3 
gauge covariant Wigner transform in the case of gauge theories) and often requires a quasiparticle approximation ^ 23j . 
The rationale behind a Wigner transform of a nonequilibrium Green's function is the assumption of a wide separation 
between the microscopic (fast) and relaxational (slow) time scales, typically justified in a weakly coupled theory. A 
recent derivation of transport equations for a hot QCD plasma along these lines has recently been reported in [ pO[ , 
however the collisional terms obtained in the quasiparticle and relaxation time approximations turn out to be infrared 
divergent. 

Thus, the importance of a fundamental understanding of transport in quantum field theory from first principles, with 
the direct application to the experimental aspects of the search for the QGP justifies the study of transport phenomena 
from many different perspectives. In this article we present a novel method to obtain quantum kinetic equations 
directly from the underlying quantum field theory implementing a dynamical renormalization group resummation. 
Such approach has been recently introduced to study the relaxation of mean-fields of hard charged scalars in a gauge 
theory . This method allowed to obtain directly in Ref. the anomalous relaxation of hard charged excitations 
in an Abelian gauge theory ||2^ ], providing an interpretation of infrared divergent damping rates ]2^ ] in terms of 
non-exponential relaxation and pointed to a shortcoming in the interpretation of quasiparticle relaxation in terms of 
complex poles in the propagator. Infrared divergences associated with the emission and absorption of long-wavelength 
gauge bosons are ubiquitous in gauge theories. Thus, this novel approach is particularly suitable to study transport 
phenomena in gauge theories. 

Goals and Strategy: The goals of this article are to provide a novel and alternative derivation of quantum kinetic 
equations directly from the microscopic quantum field theory in real time and apply this program to several relevant 
cases of interest. We consider scalar theories describing pions and sigma mesons and gauge theories. This approach 
allows to include consistently medium effects, such as nonequilibrium generalizations of the hard thermal loop resum- 
mation, describes anomalous relaxation and reveals the proper time scales for relaxation directly in real time. There 
are several advantages that this program offers as compared to other approaches to transport phenomena: 

• (i) It allows to study the crossover between different relaxational behavior in real time. This is relevant in the 
case of resonances where the medium may enhance threshold effects. 

• (ii) It describes non-exponential relaxation in a clear manner and treats threshold effects consistently, providing 
a real-time interpretation of infrared divergent damping rates in gauge theories, 

• (iii) It provides a systematic field-theoretical method to include higher order corrections and allows to incorporate 
self-consistently medium effects such as for example a resummation of hard thermal loops [p7[-^9[ that are 
necessary to determine the relevant degrees of freedom and their microscopic time scales. 

• (iv) It resolves the issue of pinch singularities that often appear in calculations of physical quantities out of 
equilibrium. 

The strategy to be followed is a generalization of the methods introduced in Ref. |Q but adapted to the description of 
quantum kinetics. The starting point is the identification of the distribution function of the quasiparticles which could 
require a resummation of medium effects (the equivalent of hard thermal loops pTj-pOf). The equation of motion for 
this distribution function is solved in a perturbative expansion in terms of nonequilibrium Feynman diagrams. The 
perturbative solution in real time displays secular terms, i.e., terms that grow in time and invalidate the perturbative 



2 



expansion beyond a particular time scale (recognized a posteriori to be the relaxational time scale). The dynamical 
renormalization group implements a systematic resummation of these secular terms and the resulting renormalization 
group equation is the quantum kinetic equation. 

The validity of this approach hinges upon the basic assumption of a wide separation between the microscopic and 
the relaxational time scales. Such an assumption underlies every approach to a kinetic description and is generally 
justified in weakly coupled theories. Unlike other approaches in terms of a truncation of the equations of motion for 
the Wigner distribution function, the main ingredient in the approach presented here is a perturbative diagrammatic 
evaluation of the time evolution of the proper distribution function in real time improved via a renormalization 
group resummation of the secular divergences. 

An important bonus of this approach is that it illuminates the origin and provides a natural resolution of pinch 
singularities | |3C| , |3l| found in perturbation theory out of equilibrium. The perturbative real-time approach combined 
with the renormalization group resummation reveal clearly that these are indicative of the nonequilibrium evolution 
of the distribution functions. In this framework, pinch singularities are the manifestation of secular terms. 

The article is organized as follows: In Sec. II we summarize the main ingredients of nonequilibrium field theory 
to establish the perturbative framework. In Sec. Ill we study the familiar case of a scalar field theory, including 
in addition the nonequilibrium resummation akin to the hard thermal loops to account for the effective masses in 
the medium and therefore the relevant microscopic time scales. In Sec. IV we discuss in detail the main features of 
the dynamical renormalization group approach to quantum kinetics, compare it to the more familiar renormalization 
group of Euclidean quantum field theory and provide an easy-to-follow recipe to obtain quantum kinetic equations. 
In Sec. V we apply these techniques to obtain the kinetic equations for cool pious and sigma mesons in the 0(4) 
linear sigma model in the chiral limit. In relaxation time approximation we obtain the relaxation rates for pious and 
sigma mesons. This case allows us to highlight the power of this approach to study threshold effects on the relaxation 
of resonances, in particular the crossover between two different relaxational regimes as a function of the momentum 
of the resonance. This aspect becomes phenomenological important in view of the recent studies by Hatsuda and 
collaborators |Q that reveal a dropping of the sigma mass near the chiral phase transition and an enhancement of 
threshold effects with potential observational consequences in heavy ion collisions. 

In Sec. VI we study the relaxation of charged quasiparticles in the full range of momenta in SQED. This theory 
has the same hard thermal loop structure at lowest order as QED and QCD p3|-^ and shares many features of 
these theories such as the lack of magnetic screening mass. In particular, in this Abelian case we provide a gauge 
invariant description of the quasiparticle distribution function, thus bypassing the complications associated with the 
gauge covariant Wigner transforms of the charged field Green's function. The hard thermal loop resummation p7|-|29|] 
is included in the scalar as well as in the gauge boson spectral densities. We find that the exchange of HTL resummed 
longitudinal photons leads to exponential relaxation but the exchange of dynamically screened transverse photons 
leads to anomalous relaxation, thus leading to a crossover behavior in the relaxation of the distribution function as a 
function of the momentum of the charged particle. The real-time description of relaxation advocated in this article 
bypasses the ambiguities associated with an infrared divergent damping rate 1^,0. In Sec. VII we discuss the issue 
of pinch singularities found in calculations in nonequilibrium field theory and establish the equivalence between these 
and secular terms in the perturbative expansion, these singularities are thus resolved via the resummation provided 
by the dynamical renormalization group. 

We summarize our results and discuss further implications and future directions in the conclusions. 



II. REAL-TIME NONEQUILIBRIUM TECHNIQUES 

The field theoretical methods to describe nonequilibrium processes have been studied at length in the literature to 
which the reader should refer for a more detailed presentation p0| , |37| "^. Here we only highlight those aspects and 
details that are necessary for our purposes. 

The basic ingredient is the time evolution of density matrix prepared initially at time t = to, which leads to the 
generating functional of nonequilibrium Green's functions in terms of a path integral defined on a contour in the 
complex time plane. 

The contour has two branches running forward and backward in the real time axis corresponding to the unitary 
evolution operator forward in time that pre-multiplies the density matrix at to and the hermitian conjugate that 
post-multiplies it and determines evolution backwards in time. The initial density matrix determines the boundary 
conditions on the propagators. 

This is a standard formulation of nonequilibrium quantum field theory known as Schwinger-Keldysh or closed-time- 
path (CTP) ||3C| , |3^^ "|43|| . Fields defined on the forward and backward branches are labeled respectively with "-f" and 
"— " superscripts and are treated independently. Introducing sources on the CTP contour, one can easily construct the 
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nonequilibrium generating functional, which generates nonequilibrium Green's functions through functional derivatives 
with respect to sources much in the same manner as the usual formulation of amplitudes in terms of path integrals. 
The path integral along the CTP contour is in terms of the effective Lagrangian defined by 

/:noncq[*+, = - c[^-] , (2.1) 

where denotes the corresponding Lagrangian in usual field theory and 4* denotes any generic (bosonic or 

fermionic) field. The advantage of the path integral representation with the above nonequilibrium, effective La- 
grangian is that it is straightforward to construct diagrammatically a perturbative expansion of the nonequilibrium 
Green's functions in terms of modified nonequilibrium Feynman rules. These nonequilibrium Feynman rules are as 
follows. 

• (i) The number of vertices is doubled: Those associated with fields on the "+" branch are the usual interaction 
vertices, while those associated with fields on the "— " branch have the opposite sign. 

• (ii) There are four propagators corresponding to the possible contractions of fields among the two branches. 
Besides the usual time-ordered (Feynman) propagators which are associated with fields on the branch, 
there are anti-time ordered propagators associated with fields on the "— " branch and the Wightman functions 
associated with fields on different branches. 

• (iii) The combinatoric factors of the Feynman diagrams are the same as those in the usual calculation of S- matrix 
elements in field theory. 

For a scalar (bosonic) field $(a;), the spatial Fourier transform of the nonequilibrium propagators are defined by 
(the extension to the case of a gauge or fermionic field is straightforward) 

G>(t,i') = j J d^xe-*''-"($(x,t)$(0,i')) , (2.2a) 
G<{t,t')^i [ d^xe^'^'' UiO,t')<^>{^,t)) , (2.2b) 



G++{t, t') - G> (i, t')e{t - t') + G< {t, t')e{t' - t) , (2.2c) 

G^-{t,t')^G>{t,t')e{t' -t)+G<{t,t')e{t~t') , (2.2d) 

G+-{t,t') = G<{t,t') , (2.2e) 

G^ + (tX)^G>{tX) , (2.2f) 

where (• • •) denotes the expecta tion value with respect to the initial density matrix. From the definitions of the 
nonequilibrium propagators eqs. ( |2.2| ), it is clear that they satisfy the identity: 

G++[t,t') + G^-[t,t') - G+-{t,t') - G^ + {t,t') - . (2.3) 

The retarded and advanced propagators are defined as 

GR.k(t, t') = G++{t, t') - G+-{t, t') = [G> {t, t') - G< {t, t')] e{t - t') , 

GA,k(i,i') - G++{t,t') - G^ + {t,t') = [G<{tX) ~ G>{t,t')\ e{t' - t) , 



which are useful in the discussion of the pinch singularities discussed in a later section (see Sec. VII). 

It now remains to specify the initial state. If we were considering the situation in equilibrium the natural initial 
density matrix would describe a thermal initial state for the free particles at temperature T. The density matrix 
of this initial state is p = exp^—Ho/T)^ where Hq is the free Hamiltonian of the system, and the time evolution 
is with the full interacting Hamiltonian. This is tantamount to switching on the interaction at t = to- If the full 
Hamiltonian does not commute with Hq the density matrix evolves out of equilibrium for t > tg. This choice of 
the thermal initial state for the free particles determines the usual Kubo-Martin-Schwinger (KMS) conditions on the 
Green's functions: 

G<{t,t')^G>{t~zf3,t') . (2.4) 
Perturbative expansions are carried out with the following real-time equilibrium free quasiparticle Green's functions: 
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G>{t,t') - ^ {[1 + nBK)] e-'^"''^*-*') +nBK) e'">'(*-*')} , (2.5a) 
G<(t,t') = ^{nB(a;k)e-^->'(*-*') + [l + nBK)]e'"''(*'*')} , (2.5b) 



V/c^ + , UBioj) = [exp(/3w) - , (2.5c) 



where (here and henceforth) k = |k|, and m is the mass of the field and ^^(a;) is the equiUbriiim Bose-Einstein 
distribution function. 

In a hot and/or dense medium the definition of the quasiparticles whose distribution function we want to study 
may require a resummation scheme such as for example that of hard thermal loops generalized to noncquilibrium 
situations. In these cases, the Hamiltonian is rearranged in such a way that part of the interaction is self-consistently 
included in the part of the Hamiltonian that commutes with the quasiparticle number operator, call it for convenience 
Hq and specific counterterms are included in the interacting part Hi to avoid double counting. 

As we are interested in obtaining an equation of evolution for a quasiparticle distribution function, the most natural 
initial state corresponds to a density matrix that is diagonal in the basis of free quasiparticles, i.e. that commutes 
with Hq. This initial density matrix is then evolved in time with the full Hamiltonian, and if the interaction does not 
commute with Hq the distribution function of these quasiparticles will evolve in time. 

The distribution function nk(to) is the expectation value of the operator that counts these quasiparticles in the 
initial density matrix. Under the assumption that the initial density matrix is diagonal in the basis of this quasi- 
particle number, perturbative expansions are carried out with the following noncquilibrium free quasiparticle Green's 
functions: 

G>(t,t')= 2^{[l + ^k(io)]e-^'^''(*-*')+nk(to)e'"''(*'*'^} , (2.6a) 

G< [t, t') = ^ {nk(to) e-'^-^'-''^ + [1 + nk(to)] e'"''^*-*') } , (2.6b) 

where Wk is the dispersion relation for the free quasiparticle. In this picture the width of the quasiparticles arises from 
their interaction and is related to the relaxation rate of the distribution function in relaxation time approximation. 
This point will become more clear in the sections that follow where we implement this program in detail. 



Finally, it is easy to check that the (bosonic) free quasiparticle Green's functions, eqs. (2.6) and (2.5), satisfy 

G>{t,t')^G<{f,t), (2.7) 
which will be useful in our following calculations. 



III. SELF-INTERACTING SCALAR THEORY 

We begin our investigation with a self-interacting scalar theory. The Lagrangian density is given by 

m^\{d,^f -\<<^' ^ (3.1) 

where mo is the bare mass. 

As mentioned in the introduction, the first step towards understanding the kinetic regime is the identification of 
the microscopic time scales in the problem. In a medium, the bare particles are dressed by the interactions becoming 
quasiparticles. One is interested in describing the relaxation of these quasiparticles. Thus the important microscopic 
time scales are those associated with the quasiparticles and not the bare particles. If a kinetic equation is obtained in 
some perturbative scheme, such a scheme should be in terms of the quasiparticles which already implies a resummation 
of the perturbative expansion. This is precisely the rationale behind the resummation of the hard thermal loops in 
finite temperature field theory |p7|-p9| and also behind the self-consistent treatment 

In a scalar field theory in equilibrium such a self-consistent resummation can be implemented by writing in the 
Lagrangian 

ml = miff + 6m^ , (3.2) 

where moff is the renormalized and temperature dependent quasiparticle thermal effective mass which enters in the 
propagators, and Sm^ is a counterterm which will cancel a subset of Feynman diagrams in the perturbative expansion 
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and is considered part of the interaction Lagrangian. As shown in Ref. |4J| for the scalar field theory case, this method 
implements a resummation akin to the hard thermal loops in a gauge theory |27-2^]. Parwani showed EJl that this 
resummation is effectively implemented by solving the following self-consistent gap equation for m^ff [pUii 




($2) = 



(27r)3 



2wk 



(3.3) 



with tJk — y/k'^ + m^^jj-. The divergences (quadratic and logarithmic in terms of a spatial momentum cutoff) in the 
zero-temperature part of eq. (3.3) can be absorbed into a renormalization of the bare mass by a subtraction at some 
renormalization scale. A convenient choice corresponds to a renormalization scale at T = and m{T = 0) = m is the 
zero-temperature mass. 

For T ^ rrics, the solution of the gap equation is given by 



2 ^ 
m + - 
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An 



O 



.In 



T 



In particular, for T 3> a/A T 3> to, we can neglect the zero-temperature mass m and obtain 



'■eS 
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(3.4) 



(3.5) 



In the massless case, TOcs serves as an infrared cutoff for the loop integrals | |44| , [46| . The leading term of eq. (3.5) 
provides the correct microscopic time scale at large temperature. 

We note that this renormalized and temperature dependent mass determines the important time scales in the 
medium but is not the position of the quasiparticle pole (or, strictly speaking, resonance). 

When the temperature is much larger than the renormalized zero temperature mass, the hard thermal loop resum- 
mation is needed to incorporate the physically relevant time and length scales in the perturbative expansion. For a 
hard quasiparticle fc ~ T, while for a soft quasiparticle k < y/X T, hence the longest microscopic time scale of the 
system is imicro ~ l/VXT ~ 1/TOeff. 



A. Quantum kinetic equation 



In this subsection we obtain the evolution equations for the distribution functions of quasiparticles. For this we 
consider an initial state out of equilibrium described by a density matrix that is diagonal in the basis of the free 
quasiparticles, but with nonequilibrium distribution functions. If the medium is hot these quasiparticles will have an 
effective mass nice which will result from medium effects, much in the same manner as the temperature dependent 
thermal mass in the equilibrium situation described above. This mass will be very different from the bare mass toq in 
the absence of medium effects and must be taken into account for the correct assessment of the microscopic time scales. 
Thus, we write the Hamiltonian in terms of the in medium dressed mass TOcff and a counterterm Sm^ = TOq — m^g 
which will be treated as part of the perturbation and required to cancel the mass shifts consistently in perturbation 
theory. This is the nonequilibrium generalization of the resummation described above in the equilibrium case. We 
emphasize that m^g depends on the initial distribution of quasiparticles. This observation will become important 
later when we discuss the time evolution of the distribution functions and therefore of the effective mass. 

We write the Hamiltonian of the theory as 



H = Ha + Hi 



int 



Ho = lJ d?x [n^ + (V$)2 + mlg $2] ^ 



4! 2 



(3.6a) 
(3.6b) 

(3.6c) 



where n(x, t) = $(x, t) is the canonical momentum, and the mass counterterm has been absorbed in the interaction. 
Here and henceforth, a dot denotes derivative with respect to time. The free part of the Hamiltonian Hq describes free 
quasiparticles of renormalized finite-temperature mass TOeff and is diagonal and Gaussian in terms of free quasiparticle 
creation and annihilation of operators a^(k) and a(k). 

With this definition, the lifetime of the quasiparticles will be a consequence of interactions. In this manner, the 
nonequilibrium equivalent of the hard thermal loops (in the sense that the distribution functions are non-thermal) 
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which in this theory amount to local terms, have been absorbed in the definition of the effective mass. This guarantees 
that the microscopic time scales are explicit in the quasiparticle Hamiltonian. 

As discussed in the previous section, we consider that the initial density matrix at time t — to, is diagonal in the 
basis of free quasiparticles, but with out of equilibrium initial distribution functions nk(to)- The Heisenberg field 
operators at time t are now written as 



$(x,t) 

n(x,<) 



/(^07jn(k,i)e''^-, nik,t)^z^[aH-k,t)-a{k,t)] , 



(3.7a) 
(3.7b) 



where a^(k, t) and a{k,t) are, respectively, creation and annihilation operators at time t and tiik = y'X^ + m^. The 
expectation value of quasiparticle number operators nk{t) can be expressed in terms of the field $(k, t) and the 
conjugate momentum n(k, t) as follows 

nk(t) = {a'^{k,t)a{k,t)) 

= ^{ <n(k, t) n(-k, t)) + ujI ($(k, t) $(-k, t)) 



(<i>(k,t)n(-k,i))-(n(k,t)<i>(~k,i)) 



(3.8) 



where the bracket (• • •) means an average o ver the Gaussian density matrix defined by the initial distribution functions 
'T'k(io)- The time-dependent distribution (B.8) is interpreted as the quasiparticle distribution function. 
The interaction Hamiltonian in momentum space is given by 



^^int — 



A 1 



1 f 5m f 

^ y n '^^^ ^) ^^(qi + q2 + qs + q4) + — J d\ $(q, t) $(-q, t) 



Taking the derivative of nk{t) with respect to time and using the Heisenberg field equations, we find 



(3.9) 



1 

2wk 



[<f>3(k,t)]H(-k,t) + n(k,i) [$3(k,i)] 



+ Sm^ ( $(k, t) n(-k, t) + H(k, t) $(-k, t) 



where we use the compact notation: 

[*'(k,t)] = J d\i d\2 d^'qs $(qi,<) $(q2,t) $(q3,t) (k - qi - qa - qa) 



(3.10) 



(3.11) 



In a perturbative expansion care is needed to handle the canonical momentum [H(k, <) = $(k, <)] and the scalar 
field at the same time because of Schwinger terms. This ambiguity is avoided by noticing that 

(H(k, t) [$3(_k, t)]) = Tr [p{to) Hk(t) [<i>3(-k, t)]] 

d 



^ Im — Tr [[<i>3(-k, t)]+ p{to) (k, t')] 

= |;([<i>3(-k,t)]+<i>-(k,o) 



(3.12) 



where we used the cyclic property of the trace and the "±" superscripts for the fields refer to field insertions obtained 
as variational derivatives with respect to sources in the forward (+) time branch and backward (— ) time branch in 
the nonequilibrium generating functional. 

We now use the canonical commutation relation between H and <i> and define the mass counterterm Sm? =^ A A/6 
to write the above expression as 
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A 



+3i 



(27r)3 



2( [$3(k, (-k, <') ) + A ( $+(k, (-k, t') ) + ( $+(k, i')*- (-k, i) 



$+(q,t)a>-(-q,i), 



(3.13) 



The right-hand side of eq. ( |3.13 ) can be obtained perturbatively in weak couphng expansion in A. Such a perturbative 
expansion is i n ter ms of the nonequihbrium vertices and Green 's functions eqs. (^|^) with the basic Green's functions 
given by eqs. ( ^.6[ ). At order 0{X) the right hand side of (3.13) vanishes identically. This is a consequence of the fact 
that the initial density matrix is diagonal in the basis of free quasiparticles. 



Figs, la-lc display the contributions up to two loops to the kinetic equation (3.13). The tadpole diagrams, depicted 
in Figs, la and lb as well as the last term in eq. (3.13) are canceled by the proper choice of A. 

An important point to notice is that these Green's functions include the proper microscopic scales as the contribution 
of the hard thermal loops have been incorporated by summing the tadpole diagrams. The propagators entering in 
the calculations are the resummed propag ators . The terms with A are required to cancel the tadpoles to all orders. 

Thus, from the formidable expression ( |3.13| ) only the first term remains after A is properly chosen in order to 
cancel the tadpole diagrams. This requirement guarantees that the mass in the propagators is the effective mass that 
includes the microscopic time scales. Hence, we find that the final form of the kinetic equation is given by 



A d 
6u!k dt' 



:$3(k,i)]+$-(~k,t') 



(3.14) 



with the understanding that no ta dpole diagrams contribute to the above equations as they are automatically canceled 
by the terms containing A in eq. ( 3.13 ). 

To lowest order the condition that the tadpoles are canceled leads to the following condition on A 



A = -3 



d^q l + 2nq(to) 
(27r)3 2uj^ 



therefore the effective mass is the solution to the self-consistent gap equation 



'cff 



d^q l + 271q(^o) 

(27r)3 2lJ^ 



q- 



(3.15) 



(3.16) 



We see that the requirement that the term proportional to A in the kinetic equation cancels the tadpole contributions 
is equivalent to the hard thermal loop resummation in the equilibrium case |Q] and makes explicit that m^g is a 
functional of the initial nonequihbrium distribution functions. 

As will be discussed in detail below, such an expansion will be meaningful for times t <C irci = |"-k(i)/"-k(i)|, where 
irci is the relaxational time scale for the nonequihbrium distribution function. For small enough coupling we expect 
that irci will be large enough such that there is a wide separation between the microscopic and the relaxational time 
scales that will warrant such an approximation (see discussion below). 



To two- loop order, the time evolution of the distribution function that follows from eq.(3.14) is given by 



hk(t) 



d'qi 



dh2 



3 2wk ] (27r)32cc;q, (2^)32c^q, (27r)32cOq3 

7Vl(to) COs[(cJk + t^qi + t^q2 + ^^qa ) (* - t")] + 37V2(to) COs[(cJk + ^qi + ^q2 - ^q3)(* - *")] 



- / dt" (27r)353(k-qi-q2-q3 

53 J t o 



-f 37V3(to)cOs[(wk - Wqi -Wq^ + UJc^:,){t - t")] +7V4(to) COs[(tJk - Wqi 



"q2 



io^At-t")] 



where 



Ui{t) = [l + nk(<)][l + nq,(i)][l + nq,(i)][l + nq3(i)] -nk(t) n^,{t) n^,{t) n^,{t) , 

M2(t) = [1 + nk(<)][l + nq,(i)][l + nq,(i)] n^,{t) - n^{t) n^,{t) n^,(t) [1 + n^,{t)] , 

N:i(t) = [l + nk(<)] «qi(t) n^,{t) [I ^ n^,{t)] - n^{t) [1 + 7iq,(t)][l + nq,(<)] n^,(t) , 

Ni(t) = [1 + n^{t)] Tlqi (t) riq, (t) 71q3 (t) - rH,{t) [1 + nq, {t)] [I + Tlq, (t)] [I + U^, (t)] . 



(3.17) 



(3.18a) 
(3.18b) 
(3.18c) 
(3.18d) 
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The kinetic equation (3.17) is retarded and causal. The different contributions have a physical interpretation in terms 
of the 'gain minus loss' processes in the plasma. The first term describes the creation of four particles minus the 
destruction of four particles in the plasma, the second and fourth terms describe the creation of three particles and 
destruction of one minus destruction of three and creation of one, the third term is the scattering of two particles off 
two particles and is the usual Boltzmann term. 

Since the propagators entering in the perturbative expansion of the kinetic equation are in terms of the distribution 
functions at the initial time to, the time integration can be done straightforwardly leading to the following equation: 

^dc.7^kk;M(^o)]""t^"""•^^^'"*"^^ 



nk(i) 



where TZ[u>, k; A/i(to)] is given by 
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d^2 



(3.19) 



)32a;q, (27r)32wq, (2^)320;^ 



(27r)3j3(k-qi-q2-q3) 



S{U! + LUq^ + + Wqg ) Afl{tQ) + 3 d{uj + UJq^ + UJ^^ - Wq., ) 7V2 {to) 



+ 3S{UJ - OJqi - CJq^ + UJq^) Af^{tQ) + S{uJ - UJq^ ~ CJq^ - Wqg ) A/'4 (io) 



(3.20) 



We are now ready to solve the kinetic equation derived above. Since 7?.[a;, k; A/i (io)] is fixed at initial time to, 
eq.( 3.19| ) can be solved by direct integration over t, thus leading to 



rik(t) = nk(to) + y / dcJ7^[cJ,k;M(^o)] 



1 - cos[(a; - LJk)(t ~ to)] 
Tr{uj - Wk)2 



(3.21) 



This expression gives the time evolution of the quasiparticle distribution function to lowest order in perturbation the- 
ory, but only for e arly times. To make this statement more precise consider the limit t ^ to m the expression between 
brackets in (^.2l|) which can be recognized from Fermi's Golden Rule of elementary time-dependent perturbation 
theory 



lim 

t — tn — *oo 



1 - C0s[(lJ - UJk.){t - to)] 



= [t- to) S{uj - cjk) 



Tr{uj - Wk)^ 

A more detailed evaluation of the long time limit is obtained by using the following expression 



% (1 - cosyt) p{y) TT t p(0) 
a V 



^ ^[p{v)-pm+o(\ 

a y \^ 



(3.22) 



(3.23) 



where a is a fixed positive number, p{y) is a smooth function for —a < y < oo and is regular at y = 0. Thus, provided 
that TZ[(jJ, k; A/i(to)] is finite ai lu — Wk, we find nk{t) is given by 



"-k(0 = "-k(io) 



A2 



7?.[wk, k; A/i(io)] {t — to) + non-secular terms 



(3.24) 



The term that grows linearly with time is a secular term, and by non-secular terms in (3.24) we refer to terms 
that are bound at all times. The approximation above, replacing the oscillatory terms with resonant denominators 
by t6{uj — LUk) is the same as that invoked in ordinary time-dependent perturbation theory leading to Fermi's Golden 
Rule. 

Clearly, the presence of secular terms in time restricts the validity of the perturbative expansion to a time interval 
t — to <C trci with 



trcl(k) 



3nk(to) 



A2 7^K,k;A/:,(^o)] 



(3.25) 



Since the time scales in the integral in eq.(3.21) are of the order of or shorter than tmicro ~ l/woff the asymptotic form 
given by ( 3.24 ) is valid for t — to ^ tmicm- Therefore for weak coupling there is a regime of intermediate asymptotics 
in time 



t. 



ro < t - to < irci(k) 



(3.26) 



such that (i) the corrections to the distribution function is dominated by the secular term, and (ii) perturbation theory 
is valid. 

We note two important features of this analysis: 
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• (i) In the intermediate asymptotic regime (3.26) the only explicit dependence on the initial time is in the 
secular term, since 7?,[wk, k; A/i(io)] depends on to only implicitly through the initial distribution functions. 
These observations will become important for the analysis that follows below. 



• (ii) 7?.[cjk, k; A/i(to)] given by eq. ( 3.20| ) with ( 3.18a )-( 3.18d ) evaluated at to vanishes if the initial distribution 
functions are the equilibrium ones as a result of the on-shell delta functions and the equilibrium relation 1 + 
n-Bi^q) = exp(/3wq)nB((xiq), in this case there are no secular terms in the perturbative expansion. 

To highlight the significance of the second point above in a manner that will allow us to establish contact with the 
issue of pinch singularities in a later section, we note that the secular term in eq. ( 3.24 ) corresponds to the net change 
of quasiparticles distribution function in the time interval t — to- To see this more explicitly, let us rewrite 



-i 

— 7^K,k;M(^o)] = TT- 



[1 + rikito)] T,< (wk, k; to) - nk(to) S>(wk, k; to) 



(3.27) 



where T,^ {uJu_,\i;to) — S^(ci;k, k; to) = 2 z Iml]i^(u;k, k; to) with Sjj(a;k, k; to) being the the on-shell retarded scalar 
self-energy ^ calcul ated to two-loop order in terms of the initial distribution functions nk(to). Indeed, the first and 
the second terms in ( 3.27 ), respectively, correspond to the 'gain' and the 'loss' parts in the usual Boltzmann collision 
term. Hence one can easily recognize that A^7?.[wk; k; A/i(to)]/3 is the net production rate of quasiparticles per unit 
tim^ Moreover, the absence of secular term for a system in thermal equilibrium [for which 7ik(to) — nBioJu)] is a 
consequence of the KMS condition for the self-energy in thermal equilibrium: 



S>(wk,k) = e'5"''I]<(wk,k) 



(3.28) 



B. Dynamical renormalization group: resummation of secular terms 

The dynamical renormalization group is a systematic generalization of multiple scale analysis and sums the secular 
terms, thus improving the perturbative expansion |^,^|. It was originally introduced to improve the asymptotic 
behavior of solutions of differential equations |^ , ^ to study pattern formation in condensed matter systems and has 
since been adapted to studying the nonequilibrium evolution of mean-fields in quantum field theory p9| and the time 
evolution of quantum systems ||5C| ]. 

For discussions of the dynamical renormalization group in other contexts, including applications to problems in 
quantum mechanics and quantum field theory, see Refs. l47|-|50|]. 

In this section we implement the dynamical renormalization group resummation of secular divergences to improve 
the perturbative expansion following the formulation presented in Ref. pif . 

This is achieved by introducing the renormalized initial distribution functions Up (r) , which are related to the bare 
initial distribution function np(to) via a renormalization constant Zp(r, to) by 



npito) = Zp(T,to) np(T) 



Zp(T,to) = 1 + — Z«(T,to) + 



where r is an arbi trary renormalization scale and Zp^''(T, to) will be chosen to cancel the secular term at a time scale 

(3.30) 



(3.29) 



T. Substitute eq. ( ^.29| ) into eq. ( |3.24| ), to ^(A^) we obtain 

A2 



nk(i) = rik(T) + y {z^''(r,to)nk(r)-f (t-to)7^K,k;A/■(T)] } + OiX^) 
To this order, the choice 



4'^(^'*o) = -(r-to) 7^[tJk,k;M(r)]/r^k(T) 



(3.31) 



leads to 



See Sec. 4.4 in Ref. M, especially pp. 83-84. 
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nk(0 = nk(r) + — {t-T) 7^K, k;A/:,(T)] + ©(A^) 



(3.32) 



Whereas the original perturbative solution was only valid for times such that the contribution from the secul ar te rm 
remains very small compared to the initial distribution function at time Iq, the renormalized solution eq. (3.32) is 
valid for time intervals t — t such that the secular term remains small, thus by choosing r arbitrarily close to t we 
have improved the perturbative expansion. 

To find the dependence of nk(''") on t, we make use of the fact that nk(t) does not depend on the arbitrary scale t: 
a change in the renormalization point t is compensated by a change in the renormalized distribution function. This 
leads to the dynamical renormalization group equation to lowest order 



d \^ 

^«k(T) - y 7^K,k;M(T)] = 



(3.33) 



This renormalization of the distribution functi on a lso affects the effective mass of the quasiparticles since m^jj 
is determined from the self-consistent equation ( |3.1(j ) which in turn is a consequence of the tadpole cancelation 
consistently in perturbation theory. Since the effective mass is a functional of the distribution function it will be 
renormalized consistently. This is physically correct since the in-medium effective masses will change under the time 
evolution of the distribution functions. 

Choosing the arbitrary scale r to coincide with the time t in eq. ( [3.33 ) , we obtain the resummed kinetic equation: 



n^{t) 



A2 



d\i 



3 2^k J (27r)32c.q, (2^)32c.q, (27r)32^q3 ^^^^ (k - Qi - q2 - Qa) 

X (5(wk + t^qi + Wqa + ^ci3)J^l{i) + 3 ^(Wk + CJqi + LO^^ - UJ^^)N2{t) 
+ 3 (5(wk - LOc^^ - t^qa + ^q3)-^3(i) +5{lO^- Wq^ - tjq^ - LO^^)Ni{t) , 



(3.34) 



where the Ni{t) are given in eqs.( 3.18a )-( 3.18d ). To avoid cluttering of notation in the above expression we have not 
made explicit the fact that the frequencies Wq depend on time through the time dependence of moff which is in turn 
determined by the time dependence of the distribution function. Indeed, the renormalization group resummation 
leads at once to the conclusion that the cancelation of tadpole terms by a proper choice of A requires that at every 
time t the effective mass is the solution of the time- dependent gap equation 



d^q 1 + 



(2^)3 2a;q(i) 



(3.35) 



where nq(t) is the solution of the kinetic eq.(3.34). Thus, the quantum kinetic equation tha t incl udes a nonequilibrium 
generalization of the hard thermal loop resummation in this scalar theory is given by ( [3.34 ) with the frequencies 
{t) given as self-consistent solutions of the time-dependent gap equation ( |3.35 ) and of the kinetic equation 



(3.34) 



The quantum kinetic equation (3.34) is therefore more general than the familiar Boltzmann equation for a scalar 
field theory in that it includes the proper in medium modifications of the quasiparticle masses. This approach provides 
an alternative derivation of the self-consistent method proposed in Ref. [Q . 

It is now evident that the dynamical renormalization group systematically resums the secular terms and the corre- 
sponding dynamical renormalization group equation extracts the slow evolution of the nonequilibrium system. 

For sm all de partures from equilibrium the time scales for relaxation can be obtained by linearizing the kinetic 
equation ( [3.34| ) around the equilibrium solution ai t = t^. This is the relaxation time approximation which assumes 
that the distribution function for a fixed mode of momentum k is perturbed slightly off equilibrium such that nk(to) = 
'T-B(wk) + ^'T-k(^o), while all the other modes remain in equilibrium, i.e., nk+q(io) = '^_B(wk-i-q) fo r q ^ 0. 

Recognizing that only the on-shell delt a fun ction that multiplies the scattering term M'i{t) in (3.34) is fulfilled, we 
find that the linearized kinetic equation (3.34) reads 



<5nk(i) = -7(k)(5nk(t) , 



(3.36) 



where 7(k) is the scalar relaxation rate 

A^TT f d\i d3 



7(k) 



92 



d3 



93 



(27r)353(k - qi - q2 - C\?,)5{uJl^ - t^q^ - CJq^ + tJq.,) 



2wk7 (27r)32wqi (27r)32tjq2 (27r)32tjq3 

X { [1 + UB (t^qi )] [1 + riB (Wq2 )] riB (Wq3 ) - Tlfi (tjq, ) Ub (cJq^ ) [1 + ?^i3 (^^qa )] } 



(3.37) 
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Solving eq. (3.36) with the initial condition 5n\s^(t = to) = dn\!_{to), we find that the quasiparticle distribution function 
in the linearized approximation evolves in time in the following manner 

Srikit) = (5nk(io) e'''''-^^^'-'"^ . (3.38) 

The linearized approximation gives the time scales for relaxation for situations close to equilibrium. In the case of 
soft momentum (T ^ mcfr » k) and high temperature AT^ » we obtain [|j 

trci(fc « 0) = [7(fc « 0)]-i « . (3.39) 

For very weak coupling (as we have assumed), the relaxational time scale is much larger that the microscopic one 
imicro ~ 1/moff « 1/VXT, since 

^ » 1 . (3.40) 



^micro A 

This verifies the assumption of separation of microscopic and relaxation scales in the weak coupling limit 



IV. COMPARISON TO THE USUAL RENORMALIZATION GROUP AND GENERAL STRATEGY 

In order to relate this approach to obtain kinetic equations using a dynamical renormalization group to more 
familiar situations we now discuss two simple cases in which the same type of method leads to a resummation of the 
perturbative series in the same manner: the first is the simple case of a weakly damped harmonic oscillator with a 
small damping coefficient and the second, closer to the usual renormalization group ideas is the scattering amplitude 
in a four dimensional scalar theory. 



A. The weakly damped harmonic oscillator 

Consider the equation of motion for a weakly damped harmonic oscillator: 

y + y = -ey , e < i ■ 

Attempting to solve this equation in a perturbative expansion in e leads to the lowest order solution 



yit) = A e'* 



e 



1- -t 
2 



+ c.c. + non-secular terms 



where the term that grows in time, i.e., the linear secular term leads to the breakdown of the perturbative expansion at 
time scales tbrcak oc 1/e. The dynamical renormalization group introduces a renormalization of the complex amplitude 
at a time scale r in the form A = Z{t) A{t) with Z{t) = 1 + zi(r) e + • • •. Choosing zi to cancel the secular term at 
this time scale leads to 

y{t) =. A{t) [l - I (i - r)] + c.c. . 

The solution y{t) cannot depend on the arbitrary scale at which the secular term (divergence) has been subtracted, 
and this independence dy{t)/dT = leads to the following renormalization group equation to lowest order in e 

Now choosing t — t, the renormalization group improved solution is given by 

y{t) = e-^*/2 ^^(0) e'* + c.c] . 

This is obviously the correct solution to 0(e). The interpretation of the renormalization group resummation is 
very clear in this simple example: the perturbative expansion is carried out to a time scale r ^ 1/e within which 
perturbation theory is valid. The correction is recognized as a change in the amplitude, so at this time scale the 
correction is absorbed in a renormalization of the amplitude and the perturbative expansion is carried out to a longer 
time but in terms of the amplitude at the renormalization scale. The dynamical renormalization group equation is the 
differential form of this procedure of evolving in time, absorbing the corrections into the amplitude (and phases) and 
continuing the evolution in terms of the renormalized amplitudes and phases. As we will see with the next example 
this is akin to the renormalization group in field theory. 
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B. Scattering amplitude in scalar field theory 



Consider the scalar field theory described by the Lagrangian density (3.1) defined as a field theory in four dimensions 
with an upper momentum cutoff A and consider for simplicity the massless case. The IPX four point function (two 
particles to two particles scattering amplitude) at the off-shell symmetric point is given to one-loop at zero temperature 
in Euclidean space by 

r^^Hp,p,p,p) = Xo-hl\n(-) +0{\l) , (4.1) 



2 ^' \p 

where Aq is the bare coupling, p is the Euclidean four-momentum. Clearly perturbation theory breaks down for 
A/p > e^/^o. 

Let us introduce the renormalized coupling constant at a scale k as usual as 

Ao - Zxin) X{k) , Zxin) - 1 + zi(k) X{k) + 0{X^) , 

and choose zi{k) to cancel the logarithmic divergence at a arbitrary renormahzation scale k. Then in terms of A(k) 
the scattering amplitude becomes 

T^^\p,p,p,p) = X{k) + h'{K)\n^ + 0{X') , (4.2) 

with r(^)(K, K, K, k) — A(k). The scattering amplitude does not depend on the arbitrary renormahzation scale k and 
this independence implies a dV'^^^p, p, p, p) / dn — 0, which to lowest order leads to the renormalization group equation 

^dX^i^, , (4.3) 

2 

where fix = |A^(k) -I- 0{X^) is recognized as the renormalization group beta function. Solving this renormalization 
group equation with an initial condition A (p) = A that determines the scattering amplitude at some value of the 
momentum and choosing k = p 'm eq. (4^), one obtains the renormalization group improved scattering amplitude (at 
an off-shell point) 

rW(p,p,p,p;p,A) = A(p) , (4.4) 
with X{p) the solution of the renormalization group equation 

A 



1 - (3A/2) \n{p/p) 



The connection between the renormalization group in momentum space and the dynamical renormalization group 
in real time (resumniation of secular terms) used in previous sections is immediate through the identification 

to ln(A/p) , t ^ ln{p/p) , T ^ ln(K/p) , 



which when replaced into (4.1) illuminates the equivalence with secular terms. 

This simple analysis highlights how the dynamical renormalization group does precisely the same in the real-time 
formulation of kinetics as the renormalization group in Euclidean or zero temperature field theory. Much in the same 



manner that the renormalization group improved scattering amplitude (4.4) is a resummation of the perturbative 



expansion, the kinetic equations obtained from the dynamical renormalization group improvement represent a re- 



summation of the perturbative expansion. The lowest order renormalization group equation (4.3) resums the leading 
logarithms, while the lowest order dynamical renormalization group equation resums the leading secular terms. 

We can establish a closer relationship to the usual renormalization program of field theory in its momentum shell 
version with the following alternative interpretation of the secular terms and their resummation |24| . 

The initial distribution at a time tQ is evolved in time perturbatively up to a time scale t^ + At such that the 
perturbative expansion is still valid, i.e., trei ^ At with tiei the relaxational time scale. Secular terms begin to 
dominate the perturbative expansion at a time scale At ^ tmicro with tmicro the microscopic time scale. Thus if there 
is a separation of time scales such that tj-oi 3> At ^ tmicroj then in this intermediate asymptotic regime perturbation 
theory is reliable but secular terms appear and can be isolated. A renormalization of the distribution function absorbs 
the contribution from the secular terms. The "renormalized" distribution function is used as an initial condition at 
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to + At to iterate forward in time to to + 2At using perturbation theory but with the propagators in terms of the 
distribution function at the time scale to + At. This procedure can be carried out "infinitesimally" (in the sense 
compared with the relaxational time scale) and the differential equation that describes the changes of the distribution 
function under the intermediate asymptotic time evolution is the dynamical renormalization group equation. 

This has an obvious similarity to the renormalization in terms of integrating in momentum shells, the result of 
integrating out degrees of freedom in a momentum shell are absorbed in a renormalization of the couplings and an 
effective theory at a lower scale but in terms of the effective couplings. This procedure is carried out infinitesimally 
and the differential equation that describes the changes of the couplings under the integration of degrees of freedom in 
these momentum shells is the renormalization group equation. For other examples of the dynamical renormalization 
group and its relation to the Euclidean renormalization program see Ref. [ p4[ . 

An important aspect of this procedure of evolving in time and "resetting" the distribution functions, is that in this 
process it is implicitly assumed that the density matrix is diagonal in the basis of free quasiparticles. Clearly if at 
the initial time the density matrix was diagonal in this basis, because the interaction Hamiltonian does not commute 
with the density matrix, off-diagonal density matrix elements will be generated upon time evolution. In resetting the 
distribution functions and using the propagators in terms of these updated distribution functions we have neglected 
off-diagonal correlations, for example, in terms of the creation and annihilation of quasiparticles a^(k) and a(k) upon 
time evolution new correlations of the form (a(k)a(k)) and its hermitian conjugate will be generated. In neglecting 
these terms we are introducing a coarse graining B , thus several stages of coarse-graining had been introduced: (i) 
integrating in time up to an intermediate asymptotics and resumming the secular terms neglect transient phenomena, 
i.e., averages over the microscopic time scales and (ii) off diagonal matrix elements (in the basis of free quasiparticles) 
had been neglected. This coarse graining also has an equivalent in the language of Euclidean renormalization: these 
are the irrelevant couplings that are generated upon integrating out shorter scales. Keeping all of the correlations in 
the density matrix would be equivalent to a Wilsonian renormalization in which all possible couplings are included in 
the Lagrangian and all of them are maintained in the renormalization on the same footing. 



C. Quantum Kinetics 



Having provided a method to obtain kinetic equations by implementing the dynamical renormalization group 
resummation and compared this method to the improvement of asymptotic solutions of differential equations as well 
as with the more familiar renormalization group of Euclidean quantum field theory we are now in position to provide 
a simple recipe to obtain kinetic equations from the microscopic theory in the general case: 

• (1) The first step requires the proper identification of the quasiparticle degrees of freedom and their dispersion 
relations that is frequency vs. momentum which is determined from the real part of the self-energies on shell. 
The damping of these excitations will arise as a result of their interactions and will be accounted for by the 
kinetic description. Define the number operator Nk{t) that counts these quasiparticles in phase space and 
split the Hamiltonian into a part that commutes with this number operator (non-interacting) and a part that 
changes the particle number (interacting). It is important that these particles or quasiparticles are defined in 
terms of the correct microscopic time scales by including the proper frequencies in their definition. In the case 
of scalar cj)* near equilibrium at high temperature the renormalized mass is the hard thermal loop resummed, 
such would also be the case in a gauge theory in thermal equilibrium in the HTL limit. This is important to 
determine the regime of validity of the perturbative expansion within which the secular terms can be identified 
unambiguously, i.e., the intermediate asymptotics. It is here where the assumption of wide separation of time 
scales enters. Although in most circumstances the non-interacting part is simply the free field Hamiltonian (in 
terms of renormalized masses and fields) there could be other circumstances in which the non-interacting part 
is more complicated, for example in the case of collective modes. The initial density matrix is usually assumed 
to be diagonal in the basis of this number operator but with nonequilibrium distribution functions at the initial 



time. The real-time propagators arc then given by (2.6) 



(2) Use the Heisenberg equations of motion to obtain a general equation for fik{t) with nk{t) — (iVk(i)) . Perform 
a perturbative expansion of this equation to the desired order in perturbation theory, using the Feynman rules 
of real-time perturbation theory and the propagators ( |2.6D . The resulting expression is a functional of the 
distribution functions at the initial time. The only time dependence arises from the explicit time dependence 



of the free propagators (2.6). Integrate this expression in time and recognize the secular terms. 



(3) Introduce the renormalization of the distribution functions as in ( [3.2S| ) with the renormalization constant 
-E(r) expanded consistently in perturbation theory as in (3.29). Fix the coefficients z("^(t) to cancel the secular 
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terms consistently at the time scale r. Obtain the renormahzation group equation from the r independence of 
the distribution function, i.e., dnk/dr = 0. This dynamical renormahzation group equation is the quantum 
kinetic equation. 

Corollary: The similarity with the renormahzation of couplings explored in the previous section, suggests that the 
collisional terms of the quantum kinetic equation can be interpreted as beta functions of the dynamical renormahzation 
group and that the space of distribution functions can be interpreted as a coupling constant space. The dynamical 
renormahzation group trajectories determine the flow in this space, therefore fixed points of the dynamical renormal- 
ization group describe stationary solutions with given distribution functions. Thermal equilibrium distributions are 
thus fixed points of the dynamical renormahzation group. Furthermore, there can be other stationary solutions with 
non-thermal distribution functions, for example describing turbulent behavior |5l| . 

Linearizing around these fixed points corresponds to linearizing the kinetic equation and the linear eigenvalues 
are related to the relaxation rates, i.e., linearization around the fixed points of the dynamical renormahzation group 
corresponds to the relaxation time approximation. 

We now implement the program described by steps (1-3) in several relevant cases in scalar and gauge field theories. 



V. 0(4) LINEAR SIGMA MODEL: COOL PIONS AND SIGMA MESONS 

In this section we consider an 0(4) linear sigma model in the strict chiral limit, i.e., without an explicit chiral 
symmetry breaking term: 

C[7V, a] = i (d.Tvf + \ id.af -^{n^ + a'- flf , (5.1) 

where tt — (tt^, tt^, tt^^) and f-n ~ 93 MeV is the pion decay constant. At high temperature T > Tc, where Tc ^ 
0{fT^) Q is the critical temperature, the 0(4) symmetry is restored by a second order phase transition. 

In the symmetric phase pions and sigma mesons are degenerate, and the linear sigma model reduces to a self- 



interacting scalar theory, analogous to that discussed in Sec. [II. Thus, we limit our discussion here to the low 
temperature broken symmetry phase in which the temperature T f-^. Since at low temperature the 0(4) symmetry 
is spontaneously broken via the sigma meson condensate, we shift the sigma field cr(x, f) = 5'(x, t) -f v, where v is 
temperature dependent and yet to be determined. In equilibrium v is fixed by requiring that ((7(x, t)) = to all 
orders in perturbation theory for temperature T < T^. In the real-time formulation of nonequilibrium quantum field 
theory, this split must be performed on both branches of the path integral. Along the forward (-f-) and backward (— ) 
branches the sigma field (x, t) is written as 

(T±(x,i) = CT^(x,t) +i; , 

with ((T* (x, t)) — 0. The expectation value v is obtained by requiring that the expectation value of (t(x, t) vanishes in 
equilibrium to all orders in perturbation theory. Using the tadpole method |^ to one-loop order the equation that 
determines v is given by 

V [v^ -fl + (tt^) + 3(^2)] = , (5.2) 

Once the solution of this equation for v is used in the perturbative expansion up to one loop, the tadpole diagrams 
that arise from the shift in the field cancel. This feature of cancelation of tadpole diagrams that would result in an 
expectation value of a by the consistent use of the tadpole equation persists to all orders in perturbation theory. 
Furthermore, away from equilibrium, when (tt^) and (cr^) depend on time through the time dependence of the 
distribution func tion , the tadpole condition (^^) implies that v becomes implicitly time dependent. 

A solution of ( |S.2[ ) with v ^ signals broken symmetry and massless pions (in the strict chiral limit). Therefore 
once the correct expectation value v is used, the one-particle-reducible (IPR) tadpole diagrams do not contribute in 
the perturbative expansion of the kinetic equation. Up to this order the inverse pion propagator reads 



2\ 



which vanishes for vanishing energy and momentum whenever u 7^ by the tadpole condition (5.2), hence Goldstone's 
theorem is satisfied and the pions are the Goldstone bosons. The study of the relaxation of sigma mesons (resonance) 
and pions near and below the chiral phase transition is an important phenomenological aspect of low energy chiral 
phenomenology with relevance to heavy ion collisions. Furthermore, recent studies have revealed interesting features 
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associated with the dropping of the sigma mass near the chiral transition and the enhancement of threshold effects 
with potential experimental consequences |^^ . The kinetic approach described here could prove useful to further 
assess the contributions to the width of the sigma meson near the chiral phase transition, this is an important study 
on its own right and we expect to report on these issues in the near future. 

With the purpose of comparing to recent results, we now focus on the situation at low temperatures under the 
assumption that the distribution functions of sigma mesons and pions are not too far from equilibrium, i.e., cool pions 
and sigma mesons. At low temperatures the relaxation of pions and sigma mesons will be dominated by the one-loop 
contributions, and the scattering contributions will be subleading. The scattering contributions are of the same form 
as those discussed for the scalar theory and involve at least two distribution functions and are subdominant in the 
low temperature limit as compared to the one loop contributions described below. 

Since the linear sigma model is renormalizable and we focus on finite-temperature effects, we ignore the zero- 
temperature ultraviolet divergences which can be absorbed into a rcnormalization of /„. For a small departure from 
thermal equilibrium, we can approximate (tt^) and (ct^) by their equilibrium values: 

2\ _ /-2\ _ f '^^l "s(t^q^ 



where Wq = \J + m^. The sigma mass = 2\v^ is to be determined self-consistently. In the low temperature 
limit T < we find = [1 - 0{T^ / fl)] and m„ = V2A/^[1 - ©(TV/^)]. Thus in the case of a cool linear 
sigma model where T <ti f^r, we can approximate v and by /^r and \/2A/7r, respectively. 

The main reason behind this analysis is to display the microscopic time scales for the mesons: tmicro.o- ^/'m-a- 
and tmicro,7r = with k being the momentum of the pion. The validity of a kinetic description will hinge upon the 
relaxation time scales being much longer than these microscopic scales. 

Finally, the Lagrangian for cool linear sigma model reads to lowest order 

£[7r, a] = i {d.Tvf + 1 {d,af - im>2 _ a/, {an' + ' ^ + ' (5-4) 

where we have omitted the bar over the shifted sigma field for simplicity of notation. 

Our goal in this section is to derive the kinetic equations describing pion and sigma meson relaxation to lowest 
order. The unbroken 0(3) isospin symmetry ensures that all the pions have the same relaxation rate, and sigma 
meson relaxation rate is proportional to the number of pion species. Hence for notational simplicity the pion index 
will be suppressed. We now study the kinetic equations for the pion and sigma meson distribution functions. 



A. Relaxation of cool pions 



Without loss of generality, in what follows we discuss the relaxation for one isospin component, say tt^, but we 
suppress the indices for simplicity of notation. As before, we consider the case in which at an initial time t = tg, the 
density matrix is diagonal in the basis of free quasiparticles, but with out of equilibrium initial distribution functions 
"-k(^o) and n^{to)- The field operators and the corresponding canonical momenta in the Heisenberg picture can be 
written as 



where 



P^(x, t) 



7r(k,t) 
a(k,t) 



(27r)3/2 



2k 



7r(k,t) \ 

/V(k, i) 1 ik x 

a(k,t) f ' ' 
P<x(k,t) J 


(5.5) 




(5.6a) 


-4(-k,t)] , 


(5.6b) 


4(-k,i)] , 


(5.6c) 


-4(-k,i)] , 


(5.6d) 
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with CiJk = ^JW- 
P7r(k, t) as 



T,^. The expectation value of pion number operator can be expressed in terms of 7r(k, and 
<(<) = (at(k,<)a,(k,t)) 

+ ifc[(7r(k,t)P,(-k,<)) - (P^(k,i)7r(-k,i)) 

Use the Heisenberg equations of motion, to leading order in A, we obtain (no tadpole diagrams are included since 
these are canceled by the choice of v) 



Kit) 



k 



dh 



d 



(27r)3/2 \dt' 



a+(k-q,t)^+(q,t)7r-(-k,t'; 



(5.7) 



The expectation values can be calculated perturbatively in terms of nonequilibrium vertices and Green's functions. 
To 0{\) the right hand side of eq. (5.7) vanishes identically. Fig. 2a shows the Feynman diagrams that contribute to 
order A^. It is now straightforward to show that n1_{t) reads 



1 



df 



'[Ni{t^) cos[(fc + q + c^k+q)(t - t")] + N2(ta) cos[(fc -q~ uj],+^){t - t")] 



where 



+ A/'3(to) cos[(fc -q + t^k+q)(i - t")] + A/; (to) cos[(fc + q- t^k+q)(i - i")]} , 



N,{t) = [1 + nimi + n-(t)][l + n^+q(t)] - nl{t) n-(<) <+q(t) , 
N^it) = [1 + Kit)] nl[t) <+q(t) - nl(t) [1 + n-(t)][l + <+q(t)] 

A/'aW = [1 + <W] «qW[l + <+qW] - <W [1 + <W] <+q(0 

AA4(i) = [1 + <(i)][l + n-(t)] <+q(i) - nl{t) nlit) [1 + <+q(t)] 



(5.8a) 
(5.8b) 
(5.8c) 
(5.8d) 



The different contributions have a very natural interpretation in terms of 'gain minus loss' processes. The first term 
in brackets corresponds to the process — > cr + tt + tt minus the process tr + tt + tt — > 0, the second and third terms 
correspond to the scattering tt + cr ^ tt minus tt —t it + a , and the last term corresponds to the decay of the sigma 
meson cr — > tt + tt minus the inverse process tt + tt — > ct. 

Just as in the scalar case, since the propagators entering in the perturbative expansion of the kinetic equation are 
in terms of the distribution functions at the initial time, the time integration can be done straightforwardly leading 
to the following equation: 



K{t) = >^^ j c^^^7^,[c^,k;M(^o)] 

where 7?,7r[ci.', k; A/i(io)] is given by 



sin[(w — k){t — <o)] 
7r(a; — k) 



(5.9) 



7^^[w,k;M(^o)] - -f 



fl f d\ 1 



5{uo + q + Wk+q) 7Vi(io) + (5(w - g - Wk+q) 7V2(to) 



(27r)3 gwk+q 

+ (5(W - g + Wk+q) Ms.{to) +5{UJ + q- Wk+q) Ni{to) 



Eq.(5.9) can be solved by direct integration over t with the given initial condition at t^, thus leading to 

1 — cos[(a; — k){t — to)] 



Kit) = nlito) + X' / dw7^,kk;M(^o)] 



7r(aj — k)'^ 



(5.10) 



(5.11) 



Potential secular term arises at large times when the resonant denomin ator in (5.11 ) va nishes, i.e., w ~ fc. A detailed 
analysis reveals that 7?,Tr[ti^, k; A/i(to)] is regular at w = fc, hence using ( 3.22 ) and ( p. 23 ) we find that at intermediate 
asymptotic time k{t — to) 3> 1, the time evolution of the pion distribution function reads 



Kit) = Ki^o) + X^ TZTr[k, k;J\fi{to)] [t — to) + non-secular terms , 



(5.12) 
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where 7?.7r[fc, k; A/i(io)] does not depend on explicitly. 

At this point we would be tempted to follow the same steps as in the scalar case and introduce the dynamical 
renormalization of the pion distribution function. However, much in the same manner as the renormalization program 
in a theory with several coupling constants, in the case under consideration the tt field and the a field are coupled. 
Therefore one must renormalize all of the distribution functions on the same footing. Hence our next task is to obtain 
the kinetic equations for the sigma meson distribution functions. 



B. Relaxation of cool sigma mesons 



As before, we consider the case in which at an initial time t = to^ the density matrix is diagonal in the basis of 
free quasiparticles, but with initial out of equilibrium distribution functions n'^{to) and nj(to)- Again, for notational 
simplicity we suppress the pion isospin index. The expectation value of sigma meson number operator can be expressed 
in terms of cr(k, t) and Pa(k,t) as 

<(t) = (4(k,t)a,(k,t)) 

= ^|(^.(k, <)P.(-k, t)) + e (a(k, t)a{-k, t)) 

+ zfc[(a(k,t)P,(-k,<)) - {P,{k,t)<j{-k,t)) 

Using the Heisenberg equations of motion to leading order in A, and requiring again that the tadpole diagrams are 
canceled by the proper choice of v, we obtain 



Kit) 



3XU f <Pq 



(27r)3/2 



(7r+(k - q, t)^+(q, t)a- (-k, t')) + 3 (a+(k - q, t)a+(q, t)a- (-k, t')) 



(5.13) 



where the factor 3 accounts for three isospin components of the pion field. The expectation values can be calculated 
perturbatively in terms of nonequilibrium vertices and Green's functions. To 0{X) the right hand side of eq. (5.13) 



vanishes identically. Fig. 3a depicts the one-loop Feynman diagrams that enter in the kinetic equation for the sigma 
meson to order A^. To the same order there will be the same type of two loops diagrams as in the self-interacting scalar 
theory studied in the previous section, but in the low temperature limit the two-loop diagrams will be suppressed 
with respect to the one-loop diagrams. Furthermore, in the low temperature limit, the focus of our attention here, 
only the pion loops will be important in the relaxation of the sigma mesons. A straightforward calculation leads to 
the following expression 



Kit) 



3\'f^ f d\ 1 



dt" 



2c^k J (2^)3g|k + q| 

M^ito) cos[(t^k + g + |k + q|)(t - t")] +JV^{to) cos[(L^k + g - |k + q|)(t - t")] 



TVf (io)cOs[(wk + + UJk+q)it - t")] + ^2 (to) COs[{uJi^ + CJq - UJk+q){t 

q-UJk+q){t-t")] 



+ ^^{to) COs[(c^k -q+\k + q\){t~ t")] + Ml {to) cos[(c^k - g - |k + q|)(t - <")] 
9 

H 

+ N^{to) cos[(wk + t^k+q)(t - t")] -f A/I (to) cos[(u;k 



t") 



where 



Kit) - [1 + nimi + n-(t)][l + K+M ~ Kit) nlit) <+q(t) 
Kit) = [1 + n^mi + nl{t)]nl+^{t) - <(t) n-(t) [1 + nl+^{t)] 



Af^it) = [1 + n^t)] nlit) [1 + nl^m - <(t) [1 



'k+q 



it) , 



nut) = [1 + <(i)] nl{t) nl+Jt) - nUt) [1 + <(t)][l + K+M ' 



(5.14) 



(5.15a) 
(5.15b) 
(5.15c) 
(5.15d) 



and 



18 



(5.16a) 
(5.16b) 
(5.16c) 
(5.16d) 



Although the above expression is somewhat unwieldy, the different contributions have a very natural interpretation 
in terms of 'gain minus loss' processes. In the first brackets (i.e., the pion contribution) the first term corresponds to 
the process — > cr + tt + tt minus the process ct + tt + tt — > 0, the second and third terms correspond to the scattering 
TT — > TT + CT minus tt + cr — )■ tt, and the last term corresponds to the decay of the sigma meson ct — > tt + tt minus 
the inverse process tt + tt ^ cr. Similarly, in the second brackets (i.e., the sigma meson contribution) the first term 
corresponds to the process ^ cr + tr + cr minus the process a -\- a ~\- a ^ 0, the second and third terms correspond 
to annihilation of two sigma mesons and creation of one sigma meson minus the inverse process, and the last term 
corresponds to annihilation of a sigma meson and creation of two sigma mesons minus the inverse process. 

Since the propagators entering in the perturbative expansion of the kinetic equation are in terms of the distribution 
functions at the initial time, the time integration can be done straightforwardly leading to the following equation: 



(i) = a2 / dcJ7^,[cJ,k;M(^o)] 



sin[(^ - uj]^){t - to)] 

Tt{uJ - Wk) 



(5.17) 



where 



7^<,[w,k;M(to)] 



dPq ] 1 

Wr\ ^Ik + ql 



5{u + q+\k + qDA/T (io) + ^(^^ + g - |k + q\)Af:;{tQ) 



+ 5{uj-q+\k + q|)A/?(to) + <5(c^ - g - |k + q\)Ul{t^) 



5[ljJ + CJq + t^k+q)A/'f (to) + 5{UJ + CJq - UJ\^+ci)J^2 {^o) 



+ 5[ijj - t^q + LJk+q)AA:^ (^o) + 5{uj - - a;k+q)A/X (^0 ) 



(5.18) 



Just as before TZa[i^-,'^]Afi{tQ)\ is fixed at initial time io, eq. ( 5.17 ) can be integrated over t with the given initial 
condition at ioi thus leading to 



<(0=<(io) + A2 / dcJ7^<,[a;,k;M(^o)] 



1 — COs[(clI — LO\^){t — to 



(5.19) 



At intermediate asymptotic times ma{t — to) 3> 1, potential secular term arises when lj ^ uj\^ in eq. ( 5.19| ). We notice 
that, although 7?,o-[w, k; A/i(to)] has threshold (infrared) singularities at w = ±fc, it is regular on the sigma meson 
mass-shell. This observation will allow us to explore a crossover behavior for very large momentum later. 

Since the spectral density is regular near the resonance region uj = it^k, the behavior at intermediate asymptotic 
times is given by 



n'^{t) — n'^{to) + X^7icr[ujk, k;f\fi{to)]{t — to) + non-secular terms. 



(5.20) 



We note that the perturbative expansions for the pion and sigma meson distribution functions contain secular terms 
that grow linear ly in time, unles s the system is initially prepared in thermal equilibrium. We must now rcnormalize 
both equations (5.12) and (5.20) simultaneously, since it is a field theory with two coupled fields. 

Introduce the renormalized initial distribution functions np(T) and n^lr), which are related to the bare initial 
distribution functions n^lto) and n^lto) via respective renormalization constants Z^^r^to) and Zp{T,to) by 



nl{to)=Z;iT,to)nl{T) , Z;{T,to) = ^ + ^' ^;^'Hr,to) + 

<(io) =Z^(T,to)<(T) , z;iT,to) = i + \^z;^'Hr,to) + 



(5.21a) 
(5.21b) 



where t is an arbitrary renormalization scale at which the secular terms will be canceled. The renormalization 

constants Zp^^\T,to) and Zp'"^\T,to) ar e cho sen so as to ca ncel the secular term at the arbitrary scale r consistently 
in perturbation theory. Substitute eq. (5.21) into eq. (5.12), consistently up to 0{X^) we obtain 
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Kit) = <(t) + [zI^^\t, to) Kir) + it~ to) 7^,[fc,k;M(T)]} + ©(A^) , 
Kit) = <(t) + A2 [zZ^'\T,to) Kir) + {t- to) 7e,K,k;M(T)]} + 0{X^) . 



To this order, the choices 



io) = - (t - ^o)7^. K, k; M (t)]/< (t) 



lead to 



<W + it-r) 7^,[fc,k;M(T)]+0(A4) , 

<W + (t-r) 7^,K,k;M(T)]+0(A4) . 

The independence of nj^(t) and ^^(i) on the arbitrary rcnormahzation scale r leads to the simultaneous set of dynamical 
renormalization group equations to lowest order: 

^Kir)^>^^n4k,k-M,{T)] , 

^n£(T) = A2 7^,K,k;M(T)] . 

These equations have an obvious resemblance to a set of renormalization group equations for "couplings" and 
where the right hand sides are the corresponding beta functions. 

As before, choosing the arbitrary scale r to coincide with the time t and keeping only the terms whose delta 
functions have support on the mass shells we obtain the kinetic equations describing pion and sigma relaxation: 

Kit)-^^ J ^fj^, ^gl'+q'^^''^ + + ~ <^^^ + ■ (5-22) 

3^A2/2 f d^q ^K-g-lk + ql) 



Kit) 



(27r)3 



9|k + q| 



{[l + <(t)]n-(tX+q(t)-<W + + , (5.23) 



The processes that contribute to (5.22) are depicted in Fig. 2b and those that contribute to ( ^.23 ) are depicted in 
Fig. 3b. 



C. Relaxation time approximation 



Thermal equilibrium is a fixed point of the dynamical renormalization group equations ( ^.22 ) and (5.23), i.e., a 
stationary solution of the kinetic equations. 

A linearized kinetic equation can be obtained in relaxation time approximation, in which only the mode with 
momentum k is slightly out of equilibrium whereas all the other modes are in equilibrium: 



SKit) 
SKit) 



-7.(k)<5<(t) , 
-7.(k)J<(i) , 



where (k) and 7o- (k) are respectively the cool pion and sigma meson relax ation rates which are identified with twice 
the damping rates of the corresponding field amplitudes. Linearizing (5.22) we obtain 



7^(k) 



ttAV^ f (Pq [nsiq) ~ nsiLUk+q)] 



k 

47rfe2 



(27r) 



S{k 



^k+q) 



In 



I _ g-/3(m^/4fe+fc) 



1 



-f3m'^/4k 



(5.24) 



This is a remarkable expression because it reveals that the physical processes that contribute to cool pion relaxation 
are the decay of sigma meson a —f tt -\- tt and its inverse process tt + tt — > cr. The form of (5.24) is reminiscent of the 
Landau damping contribution to the pion self-energy and in fact a simple calculation reveals this to be correct. The 
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sigma particles present in the medium can decay into pions and this increases the number of pions, but at the same 
time pions recombine into sigma particles, and since there are more pions in the medium because they are lighter the 
loss part of the process prevails producing a non-zero relaxation rate. This is an induced phenomenon in the medium 
in the very definitive sense that the decay of the heavier sigma meson induces the decay of the pion distribution 
function, it is a non-coUisional process. 

Such relaxation of cool pions is analogous to the induced relaxation of fcrmions in a fermion-scalar plasma induced 
by the decay of a massive scalar into fermion pairs . 

For soft, cool pion mode (A; <C T <C /,r), the pion relaxation rate reads 



exp 



4fcT 



(5.25) 



The exponential suppression in the soft, cool pion relaxation rate is a consequence of the heavy sigma mass. Our 
results of the pion relaxation rate are in agreement with the pion damping rate found in Ref. |54|. These results 
(accounting for the factor 2 necessary to relate the relaxation rate to the damping rate) also agree with those reported 
recently in Ref. |]55|] wherein a related and clear analysis of pion and sigma meson damping rates was presented. 
For the relaxation rate of the sigma mesons, we find 



7.(k) 



3^A^f2 [ d\ [l + nB(<z)+nB(|k + q| 

2wk 
3AV^ 



(27r)3 

2T, 
— — m 



1 



'Zlk + ql 

g-/3(i^k+fe)/2 

g-/3(wk-fc)/2 



(5(wk-9 - |k + q[ 



The first temperature-independent term in ^ai}^) is the usual zero-temperature sigma meson decay rate 
the finite temperature factors result from the same processes that determine the pion relaxation rate, i.e., cr <-> tt 
For soft sigma meson (fc ^ T ^ /^), we obtain 



(5.26) 



, whereas 



87rma 



coth 



"m„\ 

.It/ 



It agrees with the decay rate for a sigma meson at rest found in Refs. 

On the other hand consider the theoretical high temperature and large momentum limit k ^ mg- ^ T such that 
LO]^ — k<^T. In this limit the sigma meson relaxation rate ( 5.26 ) becomes logarithmic (infrared) divergent. The reason 
for this divergence is that as was mentioned below eq. (5.19), T Za{uj, k; A/"i) has an infrared threshold singularity at 
uj — k arising from the cont ributi on pr oport ional to in eq. (5. IS). In the presence of this thres hold singularity. 



we can no longer apply eqs. (3.22) and ( |3.23 ) and instead we must study the long time limit in (5.14) more carefully. 
Understanding the infiuence of threshold behavior of the sigma meson on its relaxation could be important in view 
of the recent proposal by Hatsuda and collaborators that near the chiral phase transition the mass of the sigma 
meson drops and threshold effects become enhanced with distinct phenomenological consequences. We expect to 
report on a more detailed study of threshold effects near the critical temperature in the near future. 



D. Threshold singularities and crossover 



As mentioned above, in the discussion following eq.(5.19), Tla-['-^,'ii-;J^i{to)] in (5.19) has threshold singularities at 
UJ = ±k arising from the emission and absorption of coUinear massless pions. For k ~ m^, the point at which the 
resonant denominator in ( 5.19 ) vanishes (i.e., oj — oju) is far away from threshold and TZa[u),\i;Afi{to)] is regular at 
this point (on-shcU), hence Fermi's Golden Rule ( 3.22 ) is applicable. However in the large momentum limit, when 
— > fc the point at which the resonant denominator vanishes becomes closer to threshold and such singular point 
begins to infiuence the long-time behavior. 

That this is the case can be seen in the expression for the relaxation rate ( 5.26| ) which displays a logari thmic 
(infrared) divergence as uju ^ k. A close inspection at the terms that contribute to 7?.o-[ijJ, k; A/i(io)] in eq.(5TS) 
reveals that the threshold divergence arising as a; « Wk — originates in the term proportional to Af^ito) which 
accounts for the emission and absorption of coUinear massless pions. 

In order to understand how this threshold divergence modifies the long-time behavior, let us focus on the mode 
of sigma mesons with momentum k 3> ^ T. This situation is not relevant to the phenomenology of the cool 
pion-sigma meson system for which relevant temperatures are T <^ m^r. However studying this limiting case will yield 
to important insight into how threshold divergences invalidate the simple Fermi's Golden Rule analysis leading to 
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on-shell delta functions in the intermediate asymptotic regime. This issue will become more pressing in the case of 
gauge theories studied below. 

To present this case in the simplest and clearest manner, we will study the relaxation time approximation, by 
assuming that only one mode of sigma mesons, with momentum k, is slightly displaced from equilibrium such that 
= nBi^k) + '5'^k(^o)j whereas all other pion and sigma meson modes are in equilibrium, i.e., n'^{to) = nsiq) 
for all q and n'^{to) — 7iB([x)q) for all q 7^ k. In this approximation and keeping t he on ly term that contributes to 
7?.cr[ti;, k;7Vi(io)] for uj « cjk, [i.e., the one proportional to A/"/ (to)], we find that eq. (5.19) simplifies to 



Sn^t) - <5<(to) 



1 - / (itJ7cr(a;, k) 



1 - C0S[(W - UJu.){t - to)] 



(5.27) 



with 



7<T(w,k) = — / -— r S{uj-q-\k + q\) 



2cjk 
y2f2 r 



3A^ 

STTWk 



(27r)3 
2T 



glk + ql 



1 + fc 1 - e-'5('^-fe)/2 



(5.28) 



At intermediate asymptotic times nia[t — to) 3> 1, the region lu ujk ~ k dominates the integral and in the limit 
k ^ > r we can further approximate 



T 



+ 0{uj - k) , 



(5.29) 



where T = 2T[1 — exp(— fc/T))] « 2T. The integral over uj in eq.(5.27) can be performed when j^ri^^'k-) is given by 
the first term in eq. ( 3.2£ ) and we obtain 



, , ,,, l-cos[(u;-^k)(t-to)] 3A2/|T 
dw7<,(w,k) fy- w Tit-to,k) 



TTlW - UJk 



for rricr {t — to) ^ 1 , where 





T 




tjj\^ — k 



+ ci[(cjk- fc)(i-io)] 



3in[(u;k -k){t~ to)] \ 

{LO^-k){t-to) J 



(5.30) 



with ci(x) being the cosine integral function 



ci X 



dt 



cost 



For fixed k, !F{t — to - k) has the following limiting behaviors 

Tit - to, k) - (t - to) {In [(t - to)TeT"i] + O {{uk ~ k)^{t - to)^)] for (c^k - k){t - to) « 1 , (5.31) 
J'(t-to,k) = (t- to) I In 





T 




jln 


+ ( 


LU]^ — k 



,(L^k-fc)2(t-to)^ 



for (wk -k){t~ to) > 1 



(5.32) 



where 7 — 0.5772157. . . is the Euler-Mascheroni constant. Thus, we see that there is a crossover time scale tc « 
(tJk — fc)^^ at which the time dependence of the function T{t — to, k) changes from tint for t — to ^ tc to linear in t 
for t — to^tc- In the large momentum limit, as the sigma meson mass-shell approaches threshold, this crossover time 
scale becomes longer such that an "anomalous" (non-linear) secular term of the form t In t dominates during most of 
the time whereas the usual secular term linear in t ensues at very large times. 

We can now proceed with the dynamical renormalization group to resum the secular terms. Introducing the 
renormalization constant by 



5<(to) =Z^(r,to)5<(r) , Z^(T,to) = 1 + A2^^(^)(r, to) 



and choosing 



(T,to) 



.F(r-to,k) 



(5.33) 



(5.34) 
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to cancel the secular divergences at the time scale r, we find that dynamical renormalization group equation 



dr A'Kk? dr 

leads to the following solution in relaxation time approximation 



d<5n-(r) ^ iX^flT d^r - <o, k) ^ ^ ^^^^^^ 



Sn^{t) — (5n£(io)exp 



47rfc2 



(5.36) 



In the large momentum limit, using (|5!3l| ) and (|5j|) we find that the crossover in the form of the secular terms 



results in a crossover in the sigma meson relaxation: an "anomalous" (non-exponential) relaxation will dominate the 
relaxation during most of the time and usual exponential relaxation ensues at very large times. 

This simple exercise has revealed several important features highlighted by a consistent resummation via the dy- 
namical renormalization group: 



Threshold infrared divergences result in a breakdown of Fermi's Golden Rule. The secular terms of the pertur- 
bative expansion are no longer linear in time but include logarithmic contributions arising from these infrared 
divergences. 

The concept of the damping rate is directly tied to exponential relaxation. The infrared divergences of the 
damping rate reflect the breakdown of Fermi's Golden Rule and signal a very different relaxation from a simple 
exponential. 

Whereas the usual calculation of damping rates will lead to a divergent result arising from the infrared threshold 
divergences, the dynamical renormalization group approach recognizes that these threshold divergences result 
in secular terms that are non-linear in time as discussed above. While in relaxation time approximation linear 
secular terms lead to exponential relaxation and therefore to an unambiguous definition of the damping rate, non- 
linear secular terms lead to novel non-exponential relaxational phenomena for which the concept of a damping 
rate may not be appropriate. 

This discussion of threshold singularities and anomalous relaxation has paved the way to studying the case of gauge 
theories, wherein the emission and absorption of (transverse) photons that are only dynamically screened lead to a 
similar anomalous relaxation p4[ . 

VI. HOT SCALAR QED 

In this section we study the relaxation of the distribution function of charged scalars in hot SQED as a prelude to 
studying the more technically involved cases of hot QED and QCD Hot SQED shares many of the important 
features of hot QED and QCD in leading order in the hard thermal loop (HTL) resummation psj-poj. Furthermore, 
the infrared physics in hot QED captured in the eikonal (Bloch-Nordsieck) approximation has been reproduced 
recently via the dynamical renormalization group in hot SQED p^ , thus lending more support to the similarities of 
both theories at least in leading HTL order. However, unlike hot QED and QCD there are two simplifications [^,0 in 
this theory that allows a more clear presentation of the relevant results: (i) there are no HTL corrections to the vertex 



and (ii) the HTL resummed scalar self-energy is momentum independent |33|3J]. These features of hot SQED enable 
us to probe the relaxation of charged scalars with arbitrary momentum within a simplified setting that nevertheless 
captures important features that are relevant to QED and QCD. This study is different from those in Ref. in that 
we here include the contribution from the longitudinal, Debye screened photons and discuss in detail the crossover 
between the relaxational time scales associated with the transverse and longitudinal photons for arbitrary momentum 
of the charged scalar. Furthermore, in order to provide an unambiguous definition of the distribution function, our 
study is done directly in a gauge invariant formulation. This formulation has several advantages, in that gauge 
invariance is built in from the outset and the distribution functions are defined for gauge invariant objects. 

In the Abelian theory under consideration, it is rather straightforward to implement a gauge invariant formulation 
by projecting the Hilbert space on states annihilated by the two primary first class constraints: Gauss' law and 
vanishing canonical momentum conjugate to the temporal component of the gauge field. Gauge invariant operators 
are those that commute with both constraints and are obtained systematically, finally the Hamiltonian and Lagrangian 
can be written in terms of these gauge invariant operator s ||67| , details are presented in Appendix A. The resulting 
Lagrangian is exactly the same as that in Coulomb gauge ||67|| and is given by (see Appendix A) 
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£ = 9^$^ 9^$ - 771^$^$ + ^df^Ar ■ 9^ At - cAt • Jt - c^At • At ^^<i> 

= i[$tv$ - (V$^)$] , jo = _ $t$) ^ 

where e is the gauge coupUng, At is the transverse component of the gauge field satisfying V • At(x, i) = 0, 
$ and <i>^ are charged but gauge invariant fields, and we have traded the instantaneous Coulomb interaction for 
a gauge invariant auxiliary field which should not be confused with the time component of the gauge field. 
Since we are only interested in obtaining the relaxation behavior arising from finite temperature effects we do not 
introduce the renormalization counterterms to facilitate the study, although these can be systematically included in 
our formulation p^ . Furthermore we will consider a neutral system with vanishing chemical potential. 

Medium effects are included via the equilibrium hard thermal loop resummation, hence we will restrict our study 
to the relaxation time approximation in which only one mode of the scalar field, with momentum k is perturbed 
off equilibrium while all other scalar modes and the gauge fields will be taken to be in equilibrium. Considering 
the full nonequilibrium quantum kinetic equation will require an extrapolation of the hard thermal loop program to 
situations far from equilibrium, clearly a task beyond the scope of this article. Hence the propagators to be used in 
the calculation for the modes and fields in equilibrium will be hard thermal loop resummed. 

Since for hot SQED the leading one-loop contributions to scalar self-energy is momentum independent and ^ 
O(e^T^) iQ, the leading order HTL resummed inverse scalar propagator reads (here and henceforth, we neglect the 
zero-temperature scalar mass m) 

A~\Lj,k)^u;'-k'-ml , (6.1) 
where 777^ = eT/2 is the thermal mass of the charged scalar. The dispersion relation of scalar quasiparticles to leading 



order in HTL is given by Wk = ^yW + rrl^■ Just as in the scalar case studied in section IH, the mass 777^ is included 
in the Hamiltonian and a counterterm is considered as part of the interaction to cancel the tadpole contributions. 
In terms of the free scalar quasiparticles of mass 77is, the field operators in the Heisenberg picture are written as 

^(^^t)^Jj0j,HKt)^^-, 

where 

<jy{k,t) = [a(k,t)+6t(_k,i)] , 

^{k,t)^i^[a\-k,t)-b{k,t)\ . 

The number of positively charged scalars (which at zero chemical potential is equal to the number of negatively 
charged scalars) is then given by 

77k(t) = {a\k,t)a{k,t)) 

= ^|(7r(-k,<)^t(k^t)) +^2 (0t(_k,<)</,(k,t)) 



-7Wk 



(^t(_k,<)^t(k,t))_(^(_k,<)<^(k,t)) 



We emphasize that this number operator, is a gauge invariant quantity by construction. Using the Heisenberg 
equations of motion, to lowest order in e, we obtain 

r7k(t) = 77L,k(i) + riTMit) , 

where 7iT.k(i) and nx^k{t) correspond to the longitudinal photon (plasmon) and transverse photon contributions 
respectively: 
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e 

X 



(27r)3/2 
t^+(-k,i' 



-(k-q,i)^0 (q,t) 



nT,-k{t) = — 



^^k7 (27r)3/2 



■'+(-k, i')'/'~(k-q,tMo(q, i) 



(6.2) 



+(k,i')'/'^'"(-k-q,iMr(q>i) 

X^+(q,t)^+(k-q,i)(/.t--(-k,i') 



Here k7-(q) = k — (k • q)q, A.T{k^ t) and -4o(k, t) are the spatial Fourier transforms of the gauge fields: 



At(x, t) 



(27r)3/2 



^T(k,t)e 



^o(x,t) 



(27r)3/2 



-4o(k, i) e 



(6.3) 



(6.4) 



As usual the expectation values are computed in noncquilibrium perturbation theory in terms of the real-time prop- 
agators and vertices. A detailed study of this scalar theory has revealed that there are no HTL vertex corrections in 
SQED |33y34| ] and this facilitates the analysis of the time evolution of the distribution function for soft quasiparticles. 
Since in SQED the leading order HTL contribution to the scalar propagator is a mass shift, the real-time HTL 



effective scalar propagator is given in eqs. (2.6) in terms of the quasiparticle frequency = y^^^l-mf. When the 



internal photon lines in the Feynman diagrams for the kinetic equation are soft, an HTL resummation of these photon 
lines is required [p7|- p9|j4^ ] . It is important to note that the HTL resummed photon propagators are only valid in 
thermal equilibrium since the KMS condition that relates the advanced and retarded Green's functions has been used to 
write these in terms of the spectral density. Therefore an analysis of the kinetic equation for the distribution function 
that uses the HTL resummation for the soft degrees of freedom will be restricted to the linearized, i.e. relaxation time 
approximation. A truly noncquilibrium description of the kinetic equations for charged or gauge fields will require an 
extension of the hard thermal loop program to situations far away from equilibrium, clearly such extension is beyond 
the scope of this article. Therefore the derivation of the kinetic equation for the charged scalar fields assumes that 
the photons are in equilibrium and the distribution function of the charged scalars has been displaced slightly off 
equilibrium. Fig. 4a shows the lowest order 0(62) contribution to the kinetic equation from longitudinal photons and 
Fig. 4b shows the contributions from transverse photons. 



A. Longitudinal photon contribution 

hi this gauge invariant formulation, the longitudinal photon is associated with the auxiliary field Ao(x, t) which 
is the Lagrange multiplier associated with Gauss' law constraint. Since this is not a propagating field (no canonical 
momentum conjugate exists), proper care must be taken in obtaining the Green's functions for this field. In appendix B 
we provide the details to obtain the HTL resummed real-time Green's function for this auxiliary field. 

The HTL effective propagators of the longitudinal photons are given by (see Appendix B) 



Ql.sS.^.t') = -i j d?xe-'^-^ (^o(x,<)Ao(0,i') 
^L,q(*.^') = -* / d^xe-'^-^ (^o(0,t')^o(x,t) 

QtXii, t') = ^sit - 1') + gl^{t, t')9{t - 1') + g<^it, mt' t) , 
Qr^it, t') = -\5{t t') + g>(t, t')e{t' - 1) + g<{t, t')e{t ~ t') 



,<(» 



{t,t') 



(6.5a) 

(6.5b) 

(6.5c) 

(6.5d) 
(6.5e) 



where q = |q| and 
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The HTL spectral density pL{qo, q) is given by [ p3| , p5[ 

1 ImSi(go,q)e(g'-gg) 



PL(go,q) = 

ImEL(go,q) = 
ReS]L(go,q) = 



^ [g2 +ReSi(qo,qr + [ImSi(qo, q)]' 



+ sgn(go)^L(q)'5(go -^i(q)) , 





90 


6 


9 








2 - 


6 





qo 
q 



qo + q 



qa-q 



(6.6a) 
(6.6b) 

(6.7a) 
(6.7b) 
(6.7c) 



where <j-'L(q) is the longitudinal photon pole and Z]^(q) is the corresponding (momentum dependent) residue, which 
will not be relevant for the following discussion. 

Using the above expressions for the nonequilibrium propagators, and after some tedious but straightforward algebra, 
we find nL^\^{t) to lowest order in perturbation theory ©(e^) is given by 



dPq 



2a;k J (27r)3tJk+q J-c 



/OO f^t 
dqo PL{qo,(l) / dt" 
-OO J to 

X (ujk - t^k+q)^ M (to) cos[(a;k + t^k+q + go) (i - t")] 

+ (tJk + (^k+qf J^2{to) COs[(wk - Wk+q - <7o)(i - t" 



where 



A/'i(i) = [1 + nk(0] [1 + nk+q(t)] [1 + nsiqo)] - rik(0 f^k+q(^) neiqa) , 
A/2(i) = [1 + nk(0] nk+q{t) ns(go) - »^k(i) [1 + n^+q{t)] [1 + ni3(qo)] • 

To obtain eq. (|6.^), we have used the following properties (see also Appendix B) 

PL(-?o,q) -pL(go,q) , nsi-qQ) ^ -[l + nsiqa)] ■ 



(6.8) 



(6.9a) 
(6.9b) 



(6.10) 



The different contributions have a very natural interpretation in terms of gain minus loss processes. The first term in 
brackets corresponds to the process ^ 7^ + s + s minus the process j^ + s + s ^ 0, and the second term corresponds 
to the scattering in the medium 7^ + s — > s minus the inverse process s — s- 7^ + s, where 7^ refers to the HTL-dressed 
longitudinal photon and s, s refer to the charged quanta of the scalar field <i>. 

As mentioned above the HTL resummation of the internal photon and scalar lines assume that these degrees of 
freedom are in thermal equilibrium and that the kinetic equation is valid in relaxation time approximation which 
will be assumed henceforth. Namely, we assume that at time t — to the distribution function for a fixed mode with 
momentum k is disturbed slightly off equilibrium such that n£^k(^o) = ?^s(^k) + <5?iL.k(io), while the rest of the modes 
remain in equilibrium, i.e., ni^k+q(io) = ^^(i^k+q) for q 7^ and linearize the kinetic equation in (Srix,_k- 

Since the propagators entering in the perturbative expansion of the kinetic equation are in terms of the distribution 
functions at the initial time, the time integration can be done straightforwardly leading to a linearized equation in 
relaxation time approximation. In terms of the spectral density 



2^2 /■ d^q /-c 
PL(w,k) = / — — / 

t^k J (27r)-^Wk+q 



dqo PLiqo,<i) [1 + nBiqo) + n_B (wk+q)] 

X [(O^k + t^k+q)^ S{UJ - CJk+q - qo) - (^k - t^k+q)^ S{u! + Wk+q + ^o)] , 

we obtain the time derivative of distribution function in the form 
where a = /4tt and 



(6.11) 



(6.12) 
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rL.k(t) = / dujpL{uj,k.) . (6.13) 



Upon integrating over t with the given initial condition at to leads to the form 



J to 



(6.14) 



As a consequence of the HTL resummation, the long-range instantaneous Coulomb interaction is screened with a 
Debye screening length of ©(l/eT). This results in that there are no threshold or mass shell singularities in the 
spectral density /9l(w, k) which after HTL resummation is a regular fu nctio n of oj both at threshold and on the mass 
shell uj — tJk- Therefore the analysis leading to Fermi's Golden Rule ( ^.22| ) is valid and at intermediate asymptotic 
times ms{t — Iq) ^ 1 we find a secular term that grows linear in time: 



/ rL^k(t') dt' — {t — to) pL(ijJk,k) + non-secular term 

J to 



As before applying the dynamical renormalization group to resum the secular term, one obtains the dynamical 
renormalization group (kinetic) equation 

<JnL,k(t) = -7L(k)5ni,k(t) , (6.15) 
where 7L(k) is the scalar relaxation rate corresponding to exchange of a longitudinal photon 

2TT^a f d^q (t^k + t^k+q)^ ~ , ^ 

71, (k) = / ——3 pL(uJk - 

t^k J (2ny Wk+q 

[1 + nB(cjk+q) + ?^B(t^k - t^k+q)] • (6.16) 



X 



Note that in obtaining 7L(k) we have discarded the second term in phi^'k.i k) which vanishes due to kinematics. With 
the initial condition 5nL,\i_{t — to) — <5ni.k(to)i we find that the distribution function evolves in time as 

SuL.Mit) = SnL.M e-^-('')(*-*°) . (6.17) 

Numerically, we find that 7L(fc) is a rather smooth function of k and approaches a constant value for k ^ T. The 
numerical values of 7L(fc) for static and hard scalars are, respectively, 7L(fc « 0) ~ 0.721 aT and 7i,(fc ^T) ^ 1.10 aT 
and interpolates monotonically in this range . Our results of the scalar relaxation rate due to longitudinal photon 
contribution are in agreement with the corresponding scalar damping rate found in Ref. [Q . For further comparison 
with the transverse photon contribution, we write 7L(k) in the form 

7L(k) = aTf{k) , 0.721 < /(fc) < 1.10 , (6.18) 

with /(fe) a smooth function of k. 

B. Transverse photon contribution 

We anticipate that the transverse photon contribution will lead to infrared divergences because the transverse 
photons are only dynamically screened through Landau damping in the HTL approximation [p^- p6| , ^ . Since the 
scalar is massive (ms ~ eT) the infrared region in the internal loop momenta comes solely from soft transverse photons 
with goi 9 ~ 0. In terms of the spectral density, the HTL effective nonequilibrium transverse photon propagators 
read 

V''{ci)g^,^{t,t') = I J d\e-'^-^ (A^(x,t)A^^(0,t')) , (6.19a) 

^'■''(q)^T,q(*:^') = « j d'a;e-''i "(A^j,(0,t')A^(x,i)^ , (6.19b) 
QTX^t,t') = gl^{t,t')6{t^t')+g<^^{t,t')6{t' ^t) , (6.19c) 

^T,q t') - ^T,q(i, t'Wit' 5T,q(<: t' M ~ : (S-Wd) 

gTZ(t,t')^g<>\t,t') , (6.19e) 
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where 



and 7^*-'(q) = 5^^ — (fq^jq^ is the transverse projector. Here, the HTL spectral density priqo, q) is given by 



Inii;T(Qo,q) = 
ReST(go,q) = 



1 








qo f 


12 


q \ 




\ ql 


12 


q' 



lmST{q^,Ci)e{q^ - ql) 



- Rei;T(qo,q)] + [imi;T(go,q) 

1-4 

q^ 



+ sgn((7o) ^T(q) '5(<?o " ^T(q)) , 



90 



2 + ( 1 - 1^ ) hi 



90 + g 



90 - 9 



(6.20a) 
(6.20b) 

(6.21) 
(6.22) 
(6.23) 



where ^^(q) is the transverse photon pole and Z^ici) is the corresponding (momentum dependent) residue, which 
will not be relevant for the following discussion. The important feature of this HTL spectral density is its support 
below the light cone. That is, for q^ > q^ the imaginary part of the HTL resummed photon self-energy, ImI]j-(go, q), 
originates in the process of Landau damping p^ , p5t from scattering of quanta in the medium. 

Using the above expressions for the nonequilibrium propagators and after some tedious but straightforward algebra, 
we find that nT,k(0 to lowest order in perturbation theory O(e^) is given by 



nT,k(i) 



2e^ 



d'q k2,(q) 



dqo Priqo, n) / dt" 



X 



(27r)3 tJk+q J- 

7Vi(to) cos[(tJk + ujk+q + 9o)(i - t")] 
AAj(to) cos[(a)k - t^k+q - 9o)(i - t")] 



(6.24) 



where Mijt) and Af2{t) are the same as that in eq. (3^). To obtain eq. ( 6.24 ), we have used the following proper- 
ties 0,M 



PT(-9o,q) = -pT(9o,q) , ^^B(-9o) = -[1 + ^5(90) 



(6.25) 



The different contributions have a very natural interpretation in terms of gain minus loss processes. The first term in 
brackets corresponds to the process — s- j^ + s + s minus the process 7^ -I- s + s — > 0, and the second term corresponds 
to the scattering in the medium 7^ + s ^ s minus the inverse process s — > 7^ + s, where 7^ refers to the HTL-dressed 
transverse photon and s, s refer to the charged quanta of the scalar field <&. 

As mentioned above the HTL resummation of the internal photon and scalar lines assume that these degrees of 
freedom are in thermal equilibrium and that the kinetic equation is valid in relaxation time approximation. Hence 
we assume that at time t = to the distribution function for a fixed mode with momentum k is disturbed slightly 
off equilibrium such that nr,k(^o) = '^B('^k) + 5"-T,k(^o)j while the rest of the modes remain in equilibrium, i.e., 
nT,k+q(io) = "■^(wk+q) for q ^ and linearize the kinetic equation in i5nr_k- 

Since the propagators entering in the perturbative expansion of the kinetic equation are in terms of the distribution 
functions at the initial time, the time integration can be done straightforwardly. In terms of the spectral density 



PT{c^,k) = / -—3 / 

t^k J (27r)-^ Wk+q J- 



X [1 + nsiqa) + ^^(wk+q)] S{ll> 



dqo priqo, q) 

iJ^k+q 



90 j 



we obtain the time derivative of distribution function in the form 

SnT,ic{t) = -aVuit) 5nT,kito) , 

where 

"sin[(tj - ujk)it ~ to)] 



rT,k(i) = / duj pri(^,'k.) 



7r(w - Wk) 



-co) 



(6.26) 



(6.27) 



(6.28) 
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It is to be noted that the spectral density in eq. ( |6.26 ), up to a prefactor, is the same as that studied within the 



context of the relaxation of the amplitude of a mean-field in SQED [p4| and in the eikonal approximation P5[ in QED. 
Upon integrating over t with the given initial condition at to leads to the form 



1-a [ Tu{t')dt' 



(6.29) 



As before at intermediate asymptotic times msit — to) 3> 1 if there are no singularities arising from the spectral density 
as w — > iwk, one finds a secular term linear in time. This is a perturbative signal of pure exponential relaxation at 
large times as we have discussed thoroughly in the previous sections. However, in the case under consideration the 
spectral density has an infrared singularity |^,^ and the long-time limit must be studied carefully. 

Pot ential secular terms (growing in time) could arise in the long time limit t ^ to whenever the denominators in 
( |6.28 ) vanish , i.e., for the region of frequencies uj « ±Wk. For w « Wk we see that the argument of the delta function 



in cq. (|6.26| ) vanishes in the region of the Landau damping cut of the exchanged transverse photon Qq < and 
contributes to the infrared behavior. On the other hand, for to « — Wk the delta function in eq. ( |6.26 ) is satisfied for 



qo » — 2 [Jk; and this region gives a negligible contribution to the long time dynamics. Therefore, only the first term 



in (6.2^) (with lu — ujk) contributes in the long time limit. 



This term is dominated by the Landau damping region of the spectral density of the exchanged soft photon given 



by (6.21), since for « (jjk the argument of the delta function is go + kq cos 6 /uJk and this is the region where the 
imaginary part of the HTL photon self-energy, ImI]r(go,q)) has support. The second contribution (with a; + cjk) 
oscillates in time and is always bound and perturbatively small. 

To extract the infrared behavior of the spectral density, we focus on the infrared region of the loop momenta 



with go, g ^ eT in eq._( 6.26 ) |2J,|2^. This is the region dominated by the exchange of very soft (HTL resummed) 
transverse photons p5| , ^6| | and that dominates the long time evolution of the distribution function. For go ^ 9 ^ eT, 
the contributions to the spectral density from zero temperature and massive scalars are sub-leading, therefore the 
term 1 + nB{i^k+ci) can be neglected. The only dominant contribution is from very soft quasistatic (go ^ 0) transverse 
photons for which nsiqo) ~ T/qo- 

For g <C eT the function pT {qo , q) /go is strongly peaked at go = and is well approximated by p^,|25|] 



PT(go,q) 



' ^ -^^'''^ (6.30) 



90 < 9 



7rg^ go + 1^ 



as g ^ 0, where d — Vlq^ jixe^T"^ . The remaining delta function ^(lo — Wk+q) is satisfied in the kinematical region 

qi < q < q2, with 

qi = \k - ^/uj^~m^\ , q2 ^ k + ^JZp^^^n^^ . 
The secular terms arise in the limit w ^ tJk, in this limit gi \lo — Wkl/'^k with Uk — dio-^j dk being the group veloc ity 



of the scalar quasiparticle, and g2 2k. However the region in which the above quasistatic approximation (6.3C) is 
valid corresponds to g < eT, therefore the upper momentum cutoff g2 in the integration region for g should be the 
minimum between 2k or eT. Thus for momenta k > eT the upper limit should be taken as g2 ~ eT whereas for 
k <^ eT the upper limit is g2 = 2k. 

Hence, we find that the spectral density diverges logarithmically as w — > Wk^ 

where /ik inin(wpi, k) with Wpi eT being the plasma frequency. 

As it will be seen shortly, the external momentum dependence of the upper momentum cutoff is crucial to determine 
the relaxational time scale of hard and soft scalars. At intermediate asymptotic times ujp\{t — to) ^ 1 (recall that 
u>pi ^ rus), we find |Q 

/ ^T.kit') dt' ~ 2vuT{t — to) ln[/iki'k(i ^ ^o)] + non-secular terms , (6.31) 

J to 

where /Ik = Aikexp(7 — 1) with 7 = 0.5772157. . . being Euler-Mascheroni constant. In lowest order in perturbation 
theory, the distribution functions that enters in the loops are those at the initial time. Obviously perturbation theory 
breaks down at time scales 
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t-to 



1 



1 



2avkT lii[iikVk{t — to)] 2avkTh[i{fj.]s_/2aT) 



(6.32) 



Now we apply the dynamical renormalization group to resum the anomalous secular term, (t — t^) ln[/2kWk(i — Ui)], 
in the perturbative expansion. To achieve this purpose we introduce a renormalization constant for the distribution 
function that absorbs the secular divergences at a fixed time scale t and write 



5nT,-kitQ) ^ Z{T,to)6nT,kiT) , Z(r, io) = 1 + a zi(t, to) + 



(6.33) 



and request that the coefficients z„ cancel the secular divergences proportional to a" at a given time scale r. To 
lowest order the choice 



zi{t, to) 



TTM(t')dt' 



(6.34) 



to 



leads to the rcnormalized distribution function at time t in terms of the updated distribution function at the time 
scale T 



I- a TT,k{t')dt' 



SnT,k{t) = SriT.i^ir) 



However, the distribution function SnT^k{t) cannot depend on the arbitrary renormalization scale r, this independence 
on the renormalization scale leads to the renormalization group equation to lowest order: 



dT 



SnT.kir) +aTT,kiT) 571tm{t) = 



(6.35) 



This renormalization group equation is now clearly of the form of a kinetic equation in relaxation time approximation 
with a time dependent rate. 

Now choosing the renormalization scale to coincide with the time t in the solution of ( |6.35|) as is usually done in 
the scaling analysis of the solutions to the renormalization group equations, we find that the distribution function in 
the linearized approximation evolves in time in the following manner 



SriT.i^it) = (5nT,k(to) exp 



- f rT,kit')dt' 

Jto 



(6.36) 



with the initial condition SnT.k{t = to) 5nT.k{to). In the long time limit Wpi(t — to) 3> 1, using ( 6.31 ) we find that 
the distribution function relaxes towards equilibrium as 



SnT,i^{t) « Sn^ito) exp ~'2avkT{t - to) ln[^kWk(t - to)] 



(6.37) 



Furthermore, ( |6.37 ) reveals a time scale for the relaxation of the charged scalar distribution function due to exchange 
of transverse photons trci,T = 7T^(k), with 



7T(k) 



awkT [ln(l/a) +0(1) 
2k^/ms 



for k ^ aT 
for k S, aT 



(6.38) 



Note that the transverse photon contribution to the scalar relaxation rate vanishes at zero momentum and 7r(fc) <C 
7L(fc) for k < aT, this is very similar to the behavior of the damping rate of fermions in QCD found in Rcf. [M. 



C. Relaxational crossover in real time 



The real-time description of charged scalar relaxation discussed above allows us to study the crossover between 
exponential and anomalous relaxation. Combining the longitudinal and transverse photon contributions, we obtain 
in relaxation time approximation the following time evolution of the charged scalar distribution function 

5nk(t) »'5nk(to)exp{-[7i(k)+2awkrin[A2k«k(t-to)]](i-io)} • (6.39) 

From the expression for 7L(k) given by eq.( |6.18[ ) with f{k) w 1 we find that plain exponential relaxation holds for 
2wk ln[/ik^^k(t — to)] <C 1 and ujpi(t — to) 3> 1, whereas anomalous exponential relaxation with an exponent ~ tint 
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dominates for very long times. Hence there is a crossover in the form of relaxation for the charged scalar distribution 
function at a time scale (t — to) « tc, with 

tc « !^P(iZ^ . (6.40) 

For fc <C eT we have k <^ eT and Wk 1, hence the crossover time scale is exceedingly long and the relaxation 
of the distribution function is dominated by (HTL resummed) longitudinal photon exchange and is purely exponential 
in the asymptotic regime. On the other hand for k ^ eT then /2k ~ eT and Wk ~ C'(l) and t^ ^ ijJ~^ in which case the 
relaxation is dominated by (HTL resummed) transverse photon exchange and is anomalous with an exponent tint, 
hence faster than exponential and with a relaxational time scale tici = av\^T \n{\ / a) . 



VII. SECULAR TERMS VS. PINCH SINGULARITIES 



An important difference between the approach to noncquilibrium evolution described by quantum kinetic equations 
advocated in this work and that often presented in the literature is that we work directly in real time not taking 
Fourier transforms in time. This must be contrasted with the real-time formulation (RTF) of finite temperature 
quantum field theory in which there are also four propagators and a closed-time-path contour but the propagators 
and quantities computed therefrom are all in terms of temporal Fourier transforms. In thermal equilibrium the Fourier 
representation of these four propagators for a scalar field are given by 1 40 3C| , ^ 



G++{K) 



-[G-{K)r 
1 

K"^ — m^g + it 



+ 2T:inB{\k^mK'' ~mlg) , 



G+-{K) = [0(-fco) + nB(|fco|)]^(i^' - mlg) 
G-+iK) = [0iko) + nsilkomK^ - m^^) , 



(7.1a) 

(7.1b) 
(7.1c) 



where K — (/co,k) is the four- momentum and = kg — k'^, whereas out of equilibrium the distribution functions 
are simply replaced by non-thermal ones, i.e., nsdfcol) ^ "■k(to). 
Using the integral representation of the step function 



0{t) 
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uj + ie 



one can easily show that eqs. (7T) and the ones obtained by replacing the the rma l equi libri um distributions by the 
noncquilibrium ones are, respectively, the temporal Fourier transforms of eqs. ( ^.5| ) and (2.6). The temporal Fourier 
transforms of the free retarded and advanced propagators are obtained similarly and read 



Gr/a{K) = - 



1 



K'^ - mlg ±isgn{ko)e 



(7.2) 



Several authors have poin ted out that the calculations using the CTP formulation in terms of the standard form of 
free propagators in eqs. (7T) or those obtained by the replacement of the distribution functions by the noncquilibrium 
ones, lead to pinch singularities | |3C| , |3l| , |60| "|66t . 

In a consistent perturbative expansion both the retarded and advanced propagators contribute and pinch singular- 
ities arise from the product of these, for example for a scalar field this product is of the form 



GKiK)GAiK) 



'cff 



For finite e this expression is regular, whereas when e 



< sgn(fco)e] [K"^ - mlg - i sgn(A;o)e] 

0+ it gives rise to singular products such as [5{K'^ — 



(7.3) 



as discussed in Refs. [P0| , piy60| - |66[| . Singularities of this type are ubiquitous and are not particular to scalar theories. 

A detailed analysis of these pinch terms reveals that they do not cancel each other in perturbation theory unless 
the system is in thermal equilibrium ||3^,^,^ 65 1 . Indeed, this severe problem has cast doubt on the validity or 
usefulness of the CTP formulation to describe noncquilibrium phenomena pl| . Although this singularities have 
been found in many circumstances and analyzed and discussed in the literature often, a systematic and satisfactory 
treatment of these singularities is still lacking. In Rcf. BO] it was suggested that including an in-medium width of 
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the quasiparticles to replace the Feynman's e does provide a physically reasonable solution, however this clearly casts 
doubt on the consistency of any perturbative approach to describe even weakly out of equilibrium phenomena. 

Recently some authors have conjectured that pinch singularities in perturbation theory might be attributed to a 
misuse of Fourier transforms (for a detailed discussion see [|62|-|65|). As an illustrative and simple example of these type 
of pinch singularities, these authors discussed the elementary derivation of Fermi's Golden Rule in time-dependent 
perturbation theory in quantmn mechanics. In calculating total transition probabilities there appears the square of 
energy conserving (5-function, which arises due to taking the infinite time limit of scattering probabilities. In this 
setting, such terms are interpreted as the elapsed scattering time multiplied by the energy-conservation constraint 
rather than a pathological singularity. A close look at eq. ( [7.3| ) reveals that the pinch term is the square of on-shell 
condition for the free quasiparticle, which implies a temporal Fourier transform in the infinite time limit and of the 
same form as the square of the energy conservation constraint for the transition probability obtained in time-dependent 
perturbation theory. 

By assuming that the interaction duration time is large but finite, Niegawa and Greiner and Leupold p^ ] 
showed that for a self-interacting scalar field the pinch part of the distribution function can be regularized by the 
interaction duration time a^ 



n 



r'^t)c.(i-<o)rr , (7.4) 



where "~" denotes that only the pinch singularity contribution is included, t — to is the interaction duration time and 
rj^°' is the net gain rate of the quasiparticle distribution function per unit time 



pnct 



[l + 7ik(to)]S<(cjk,k)-nk(io)S>(cJk,k) . (7.5) 



Here S]^(LJk ; k) — I]^( a;k, k ) = 2 « IniSj^(a;k, k}with Sfl(a;k,k) being the retarded scalar self-energy on mass-shell. 



Comparing (7.4) with (3.24) and (7.5) with (3.27), we clearly see the equivalence between the linear secular terms in 



the perturbative expansion and the presence of pinch singularities in the usual CTP description. In the discussion 



following eq.(3.24) we have recognized that secular terms are not present if the system is in equilibrium, much in the 
same manner as the case of pinch singularities as discussed originally by Altherr |3l|,Q . Thus our conclusion is that 
pinch singularities are a temporal Fourier transform representation of linear secular terms. 

The dynamical renormalization group provides a systematic rcsummation of these secular terms and provides the 
consistent formulation to implement the renormalization of the distribution function suggested in Refs. [ ]63| , p4| . 

Hence we emphasize that the dynamical renormalization group advocated in this article, explains the physical 
origin of the pinch singularities in terms of secular terms and Fermi's Golden Rule, and provides a consistent and 
systematic resummation of these secular terms that lead to the quantum kinetic equation as a renormalization group 
equation that determines the time evolution of the distribution function. This result justifies in a systematic manner 
the conclusions and interpretation obtained in Ref. where a possible regularization of the pinch singularities was 
achieved by including the width of the quasiparticle obtained via the resummation of hard thermal loops. 

Furthermore we emphasize that the dynamical renormalization group is far more general in that it allows to treat 
situations where the long time evolution is modified by threshold (infrared) singularities in spectral densities, thereby 
providing a resolution of infrared singularities in damping rates and a consistent resummation scheme to extract the 
asymptotic time evolution of the distribution function. The infrared singularities in these damping rates is a reflection 
of anomalous (i.e., non-exponential) relaxation as a result of threshold effects. 

The pinch singularities signal the breakdown of perturbation theory, just as the secular terms in real time, however, 
the advantage of working directly in real time is that the time scale at which perturbation theory breaks down is 
recognized clearly from the real-time perturbative expansion and is identified directly with the relaxational time scale. 
The dynamical renormalization group justifies this identification by providing a resummation of the perturbative series 
that improves the solution beyond the intermediate asymptotics. 

The resolution of pinch singularities via the dynamical renormalization group is general. As originally pointed out 
in [^l|j6C|] the pinch singularities typically multiply expressions of the form (|7.5|) which vanish in equilibrium, j ust a s 



the linear secular terms multiply similar terms in the real-time perturbative expansion, as highlighted by eqs. (3.27). 
These terms are of the typical form gain minus loss, in equilibrium they vanish, but their non-vanishing simply 
indicates that the distribution functions are evolving in time and it is precisely this time evolution that is described 
consistently by the dynamical renormalization group. 



^See eqs. (14) and (29) of Niegawa Q and eqs. (14) and (22) of Greiner and Leupold |||]. Note that eq. (14) in Ref. 
contains a typographic error. 
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VIII. CONCLUSIONS 



In this article we have introduced a novel method to obtain quantum kinetic equations via a field theoretical and 
diagrammatic perturbative expansion improved via a dynamical renormalization group resummation in real time. 
The first step of this method is to use the microscopic equations of motion to obtain the evolution equation of the 
quasiparticle distribution function, this is the expectation value of the quasiparticle number operator in the initial 
density matrix. This evolution equation can be solved in a consistent diagrammatic perturbative expansion and one 
finds that the solution for the time evolution of the distribution function features secular terms, i.e., terms that grow 
in time. In perturbation theory the microscopic and relaxational time scales are widely separated and there is a 
regime of intermediate asymptotics within which (i) the secular terms dominate the time evolution of the distribution 
function and (ii) perturbation theory is valid. A renormalization of the distribution function absorbs the contribution 
from the secular terms at a given renormalization time scale thereby improving the perturbative expansion. The 
arbitrariness of this renormalization scale leads to the dynamical renormalization group equation, which is recognized 
as the quantum kinetic equation. Linear secular terms are recognized to lead to usual exponential relaxation (in 
relaxation time approximation), whereas non-linear secular terms lead to anomalous relaxation. The dynamical 
renormalization group provides a consistent resummation of the secular terms. There are many advantages in this 
formulation: 

• It is based on straightforward quantum field theoretical diagrammatic perturbation theory, hence it allows a 
systematic calculation to any arbitrary order. It allows to include resummations of medium effects such as 
nonequilibrium generalizations of hard thermal loop resummation in the quantum kinetic equation. This is 
worked out in detail in a scalar field theory. 

• It allows a detailed understanding of crossover between different relaxational phenomena directly in real time. 
This is important in the case of wide resonances where threshold effects may lead to non-exponential relaxation 
on some time scales, and also near phase transitions where soft excitations dominate the dynamics. 

• It describes non-exponential relaxation directly in real time whenever threshold effects are important thus 
providing a real-time interpretation of infrared divergent damping rates in gauge theories. This we consider one 
of the most valuable features of the dynamical renormalization group which makes this approach particularly 
suited to study relaxation in gauge field theories in a medium where the emission and absorption of soft gauge 
fields typically lead to threshold infrared divergences. This important feature was highlighted in this article by 
studying the quantum kinetic equation for the distribution function of charged quasiparticles in SQED. 

• This method provides a simple and natural resolution of pinch singularities found often when the distribution 
functions are non-thermal. Pinch singularities are the temporal Fourier transform manifestation of the real-time 
secular terms, and their resolution is via the resummation implemented by the dynamical renormalization group. 

We have tested this new method within the familiar setting of a scalar field theory, thus reproducing previous 
results but with these new methods and moved on to apply the dynamical renormalization group to describe the 
quantum kinetics of a cool gas of pions and sigma mesons described by the 0(4) linear sigma model in the chiral 
limit. This particular example reveals a crossover behavior in the case of hard resonances because of threshold 
singularities associated with the emission and absorption of massless pions. In relaxation time approximation we find 
a crossover between purely exponential relaxation and anomalous relaxation with an exponent of the form tint which 
is faster than exponential, the crossover scale depends on the momentum of the resonance. The regime of exponential 
relaxation (in the relaxation time approximation) is described by a relaxation rate which is simple related to the 
damping rate found recently for the same model |54||5^ , |5^ . The (faster) anomalous relaxation is a novel result and 
could be of phenomenological relevance in view of recent suggestions of novel threshold effects of the sigma resonance 
near the chiral phase transition ||3^ , this possibility is worthy of a deeper study and we are currently generalizing 
these methods to reach the critical region. 

We consider that the most important aspect of this article is the study of relaxation of charged quasiparticles in 
a gauge theory. As a prelude to studying quantum kinetics in QED and QCD in this article we studied the 
case of SQED. In equilibrium, this theory shares many important features with QED and QCD in leading order in 
the hard thermal loop resummation and is a relevant model to study kinetics and relaxation in the hot Electroweak 
theory ||5^ , |3^ . This Abelian theory allowed us to begin our study by providing a gauge invariant description of the 
distribution functions, thus bypassing potential ambiguities in the definition of gauge covariant Wigner transforms 
which is the usual approach. The hard thermal loop resummation for both longitudinal and transverse photons as 
well as for the scalar is included consistently in the derivation of the quantum kinetic equation for the charged scalar 
quasiparticles in relaxation time approximation. The real-time solution of the kinetic equation for the distribution 
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function features linear and non-linear secular terms which are resummed consistently by the dynamical renormal- 
ization group. The HTL longitudinal photons are Debye screened and do not lead to infrared divergences resulting 
in purely exponential relaxation with a well defined relaxation rate. On the other hand, transverse photons are 
only dynamically screened by Landau damping and the emission and absorption of photons at right angles leads to 
infrared threshold divergences resulting in anomalous relaxation. We studied in detail the crossover between purely 
exponential and anomalous relaxation. The crossover time scale depends on the momentum and for soft quasiparticles 
exponential relaxation dominates the dynamics for a longer period of time, whereas for hard quasiparticles anomalous 
(with an exponent of the form tint) dominates the relaxation. Recent approaches to quantum kinetics including 
HTL resummations have encountered infrared divergent relaxation rates [ pO[ , the dynamical renormalization group 
reveals very clearly that this is a manifestation of non-exponential relaxation arising from threshold infrared effects 
that results in a violation of Fermi's Golden Rule. The time scales that can be extracted both from the exponential 
and the non-exponential regimes agree with those obtained by Pisarski for QCD after self-consistently including 
a width for the quasiparticle in the calculation of the damping rate |26j. Therefore, the study of this Abelian model 
has indeed offered a novel method to study relaxation in real time which is a useful arena for QCD and QED. 

We envisage several important applications of the dynamical renormalization group method primarily to study 
transport phenomena and relaxation of collective modes in gauge theories where infrared effects are important, as well 
as to study relaxational phenomena near critical points where soft fluctuations dominate the dynamics. An important 
aspect of this method is that it does not rely on a quasiparticle approximation and allows a direct interpretation 
of infrared phenomena directly in real time. Furthermore, we have established a very close relationship between 
the usual renormalization group and the dynamical renormalization group approach to kinetics. We have proved 
that the dynamical renormalization group equation is the quantum kinetic equation, the coUisional terms are the 
equivalent of the beta functions in Euclidean renormalization group. Fixed points of the dynamical renormalization 
group are identified with stationary solutions of the kinetic equation and the exponents that determine the stability 
of the fixed points are identified with the relaxation rates in the relaxation time approximation. Furthermore we have 
suggested that in this language coarse graining is the equivalent to neglecting irrelevant couplings in the Euclidean 
renormalization program. This identification brings a new and rather different perspective to kinetics and relaxation 
that will hopefully lead to new insights. 
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APPENDIX A: GAUGE-INVARIANT FORMULATION FOR SCALAR QED 

In this appendix we summarize the gauge-invariant formulation | |67| for SQED with the Lagrangian density given 
by 

where 

Fi_ii, = d^Au — di/A^ , 
D^4> = + ieAf^cl) . 

The description in terms of gauge invariant states and operators is best achieved within the canonical formulation, 
which begins with the identification of canonical field variables and constraints. These will determine the classical 
physical phase space and, at the quantum level, the physical Hilbert space. 
The canonical momenta conjugate to the gauge and scalar fields are given by 

, 

0t + ie^>t ^ 
</) - ieA°<j) . 



TT — 
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Hence, the Hamiltonian is 



H = J d^xf^ in • n + ttV + (V(/)t + ieA<j>'') ■ (V0 - ieA<j)) + i (V x A) ^ 

There are several different manners of quantizing a gauge theory, but the one that exhibits the gauge invariant 
states and operators, originally due to Dirac, begins by recognizing the first class constraints (with mutually vanishing 
Poisson brackets between constraints). From here there are two possibilities: (i) The constraints become operators in 
the quantum theory and are imposed onto the physical states, thus defining the physical subspace of the Hilbert space 
and gauge invariant operators, (ii) Introduce a gauge, converting the first class system of constraints into a second 
class (with non-zero Poisson brackets between constraints) and introducing Dirac brackets. This second possibilit y is 
the popular way of dealing with the constraints and leads to the usual gauge-fixed path integral representation |6q ] 
in terms of Faddeev-Popov determinants and ghosts. We will instead proceed with the first possibility that leads to 
an unambiguous projection of the physical states and operators. Such a method has been previously used by James 
and Landshoff within a different context |^ . 

In Dirac's method of quantization ITG] there are two first class constraints which are: 



n 



sc 



, 



(Al) 
(A2) 



with p = —i((pTr—(j)^Tr^) being the scalar charge density, eq. ( |A2| ) is Gauss's law, which can be seen to be a constraint in 
two ways: either because it cannot be obtained as a Hamiltonian equation of motion, or because in Dirac's formalism, 
it is the secondary (first class) constraint obtained by requiring that the primary constraint eq. (Al) remain constant 
in time. Quantization is now achieved by imposing the canonical equal-time commutation relations 

[AO(x,i),HO(y,t)] ^t6^x-y) , 
[A^(x,i),ff(y,<),] ^tS^^H^-y) , 

[</)t(x,i),^t(y^t)] ^^<53(x-y) . 

In Dirac's formulation, the projection onto the gauge invariant subspace of the full Hilbert space is achieved by 
imposing the first class constraints onto the states. Physical operators are those that commute with the first class 
constraints. With the above equal-time commutation relations it is straightforward to see that the unitary operator 



Ua = exp 



H°A -I- 



(fa 



performs the local gauge transformations. Thus the first class constraints are recognized as the generators of gauge 
transformations. In particular. Gauss's law operator Q is the generator of time independent gauge transformations. 
Requiring that the physical states be annihilated by these constraints is tantamount to selecting the gauge invariant 
states. Consequently operators that commute with the first class constraints are gauge invariant. 

In the Schrodinger representation of field theory, in which the field operators are diagonal, states are represented by 
wave functionals, and the canonical momenta conjugate to the field operators are represented by Hermitian functional 
differential operators. The constraints applied onto the states become functional differential equations that the wave 
functionals must satisfy: 



Mo(x) 



*[A, 



(5A(x) 



= 



ie (?!)(x) 



(x) 



0t(x) 



(5(/.t(x) 



*[A, 



(A4a) 
(A4b) 



The first equation simply means that the wave functional does not depend on Aq, whereas the second equation means 
that the wave functional is only a functional of the combination of fields that is annihilated by the Gauss's law 
functional differential operator. It is a simple calculation to prove that the fields 



35 



$(x) = (/)(x) exp 
$t(x) = (/)''' (x) exp 



ze y d3yA(y) • VyG(y ~ x) 
le J d=^yA(y) • VyG(y - x) 



are annihilated by Gauss's law functional differential equation with G(y — x) the Coulomb Green's function 

VlG{y - = S\y - . 
Furthermore, writing the gauge field in terms of transverse and longitudinal components as 

A(x) = Al(x) + At(x) , 

where 

VxAi(x) = 0, V-At(x) = 0, 

one finds 

S ^ S 



v.- 



<5A(x) 



= V. 



(A5a) 
(A5b) 

(A6) 
(A7) 
(A8) 

(A9) 



Ther efore the tr ansv erse component A^ is also annihilated by the Gauss's law operator, and A in the exponential in 
eqs. ([A5al) and (|A5b| ) can be replaced by A^ . This analysis shows that the wave functional solutions of the functional 
differential equations that represent the constraints in the Schroedinger representation are of the form 



*[A,0,0t] = *[At,$,4>^] 



(AlO) 



The fields At, ^ and are gauge invariant as they commute with the constraints. The canonical momenta conjugate 
to $ and $t are found to be 



n(x) = 7r(x) exp 
n^(x) — 7r^(x) exp 



, J d3yAi(y).VyG(y-x) 
J d3yAi(y).VyG(y-x) 



(Alia) 
(AUb) 



The momentum n canonical to A can also be written in terms of longitudinal and transverse components 

n(x) -ni(x) + nT(x) . (A12) 

It is straightforward to check that both components are gauge invariant. In the physical subspace of gauge invariant 
wave functionals, matrix elements of V-II can be replaced by matrix elements of the charge density p — ~i{^Il—^''lP). 
Therefore in all matrix elements between gauge invariant states (or functionals) one can replace 



UL{^)^-eV^J d3yG(x - y)p(y) 
Finally in the gauge invariant subspace the Hamiltonian becomes 
H = ( d 



(A13) 



j SxI^Ut ■ Ht + n + (^V$t + ieAT$t) • (v* - ieAr*) + ^ (V x Arf 
2$t<i>| _ <fxd?yp{^)G{^-y)p{y) . 



(A14) 



Clearly the Hamiltonian is gauge invariant, and it manifestly has the global U{1) gauge symmetry under which 
$ transforms with a constant phase, 11 transforms with the opposite phase, and At and 11^ is invariant. This 
Hamiltonian is reminiscent of the Coulomb gauge Hamiltonian, but we emphasize that we have not imposed any 
gauge fixing condition. The formulation is fully gauge invariant, written in terms of operators that commute with the 
generators of gauge transformations and states that are invariant under these transformations. 



36 



To obtain a gauge invariant description in Lagrangian formalism, we switch to the path integral representation for 
field theory in which the vacuum-to-vacuum amplitude is defined as 



J VAtVUtV<^>VIIV<P^VIP exp i J d^x (u^ + n^^t + Ut ■ At) - i J dtH 



(A15) 



Note that the last term in the Hamiltonian, eq. (A14), is the instantaneous Coulomb interaction, which can be traded 
for a gauge invariant auxiliary field Aq, up to an overall factor, the vacuum-to- vacuum amplitude becomes 



where 



j VAqVAtVUtV^VUV^^'VU'' exp i J d'^x (u^ + n^$t _^ jj^ . - i J dtH 
H = J d^xl^^Ur ■ Ut + tfn + (^V$t + ieAT$^) • - ^cAt*) + ^ (V x A' 



(A16) 



TO2$t$ (VAo) - eAop 



(A17) 



Since the exponent is now quadratic in the conjugate momenta, we can complete the squares and evaluate the VUt, 
Vn and VIl^ integrals to obtain 



VAqVAtV'PV'P'' exp 



■J d^xC[A„,AT,<^,^^ 



(A18) 



up to an overall factor, where £[Ao, At, 'I'^] is the gauge invariant Lagrangian, 



£[^0, At, = c)^$t _ TO2^t$ + -df^Ar ■ 9^ At - cAt ■ Jt 



(A19) 



with Jt = i[<i>^(VT<i>) - (VT<&^)<i>] and jo = - $^4). Note that Aq satisfies an algebraic equation of motion 

V^Ao = ep. 



APPENDIX B: FULL PROPAGATORS FOR AUXILIARY FIELDS 



In this appendix we derive the full real-time CTP propagators for the auxiliary field in equilibrium. We will consider 
as an example the longitudinal photon field Aq in SQED (the extension to other cases is straightforward). Since an 
auxiliary field is non-dynamical, it satisfies an algebraic equation of motion without time derivative. As a consequence, 
the free longitudinal photon propagators are local in time and there is no mixture between fields on and "— " 
branches of the CTP contour. Namely, 

^S{t-t'), (Bla) 

^^S{t-t'), (Bib) 
9 

(^o(q,i)^+(-q,t'))^ = , (Blc) 

where (• • ■)o denotes expectation value of free fields in equilibrium. 

Now we consider the full longitudinal photon propagators. Neglecting the tadpole type term e'^A^^^^ [which yields 
local (momentum independent) contribution and higher order contribution, hence is irrelevant to the one-loop result 
that we are interested in] , we have 

-Cint = eAojo , 

where jo = — i($$t — $^<i>). Straightforward diagrammatic expansions show that the following equalities hold to all 
orders in perturbation theory 



(^+(q,t)^+(-q,t'))^ = 

(A-(q,i)A"(-q,t'))g = 
(^+(q,t)A-(-q,i')) = 
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(B2a) 
(B2b) 



where (• • •) denotes the fuU equiUbrium expectation value. It is convenient to introduce the current-current spectral 
densities pj*((?o,q) and pfigo^tl) defined by 



dqo 



p>{qo,ci)9{t~t')+p<{qo,ci)e{t' -t) 



(jo^(q,ibo"("q,^'))= / dgoP<(9o,q)e-''°(*-*') 



-igo(t-t') 



Inserting a complete set of eigenstates of the full interacting Haniiltonian, one obtains the KMS condition 

p<(go,q) = e-'3*p,>('?o,q) • 
In terms of the pf{qo,<i) the full retarded longitudinal photon propagator can be written as 

A(q, t)Ao{~q, t'))^ = (^+(q, t)A+{-<l, t')) - (^+(q, t)A^ (-q, t'] 

= J|^Po(go,q)e--(-*\ 



(B3a) 
(B3b) 

(B4) 



where 



Poiqo,^) 



9 4 

(7* 



(^,q) 



qo — iu + ie 



and the KMS condition eq. (B4) is used. Thus, we obtain 



Im/9o(<7o,q) = -TT (l - e 



(9o,q) 



Again using the KMS condition eq. (B4), we can finally write the full longitudinal photon propagator as 



^+(q,i)^+(-q,0) 



l<5(t - 1') + gl^it, t')9{t - 1') + g< ^{t, t')e{t' - 1) 



where 



^L,q(^. t') = Z I dqo Inipo(go, q) [1 + ^5(^0)] e'^^^^*-*') 



^L,q(^'^') = J dqQlnypo{qo,q) nsl'Zo) 



'iqo{t — t') 



(B5) 
(B6) 



It is easy to check that the KMS condition G^qit — iP,t') — q(^'^') holds. By the same token, we obtain 



^0 (q,<Mo (-q.^O 



5{t - 1') + Gi^it, mt' -t)+ g<jt, mt t') 



L,q\ 



The imaginary part of the full retarded longitudinal photon propagator, Impo((/o, q), can be calculated in perturbation 
theory via the tadpole method |35| . To one-loop order and in HTL limit [p5|, one finds Impoigp, o) — — 7rpi(go,q), 
where Pi((7o,q) is given in eq. (|6.7| ). Finally, using the condition ^{t,t') — Qi^q{t',t) [cf. eq. (|2.7j )] and making 
change of variable go ^ ^90 in eq. (K), we find that Im/9o(— ^Oi q) = ^Ini/Oo('?o, q)7 hence /Ol(— qojq) = ^/5L(<Zo,q)- 



38 



[1] J.W. Harris and B. MuUer, Annu. Rev. Nucl. Part. Sci. 46, 71 (1996); B. MuUer in Particle Production in Highly Excited 
Matter, edited by H.H. Gutbrod and J. Rafelski, NATO ASI series B, vol. 303 (1993); H. EIze and U. Heinz, Phys. Rep. 
183, 81 (1989); H. Meyer-Ortmanns, Rev. Mod. Phys. 68, 473 (1996); H. Satz, in Proceedings of the Large Hadron Collider 
Workshop, edited by G. Jarlskog and D. Rein (CERN, Geneva), Vol. 1, page 188; and in Particle Production in Highly 
Excited Matter, edited by H.H. Gutbrod and J. Rafelski, NATO ASI series B, Vol. 303 (1993). 

B. MuUer, The Physics of the Quark- Gluon Plasma, Lecture Notes in Physics, Vol. 225 (Springer- Verlag, Berlin, Heidelberg, 
1985); L.P. Csernai, Introduction to Relativistic Heavy Ion Collisions (John Wiley and Sons, England, 1994); C.Y. Wong, 
Introduction to High-Energy Heavy Ion Collisions (World Scientific, Singapore, 1994). 

K. Geiger, Phys. Rep. 258, 237 (1995); Phys. Rev. D54, 949 (1996); Phys. Rev. D56, 2665 (1996); Phys. Rev. D46, 4965 
(1992); Phys. Rev. D47, 133 (1993); Quark-Gluon Plasma 2, edited by R.C. Hwa (World Scientific, Singapore, 1995). 

H. -T. Elze and U. Heinz, in Quark-Gluon Plasma edited by R.C. Hwa (World Scientific, Singapore, 1990); H.-T. Elze and 
U. Heinz, Phys. Rep. 183, 81 (1989). 

S. Mrowczyriski, in Quark-Gluon Plasma Q; S. Mrowczyhski and P. Daniel ewicz, Nucl. Phv s. B342, 345 (1990); S. 
Mrowczyhski and U. Heinz, Ann. Phys. (N.Y.) 229, 1 (1994); S. Mrowczyhski, |hep-ph/9805435| . 
S.A. Bass et. al. Prog. Part. Nucl. Phys. 41 225 (1998). 

I. D. Lawrie and D.B. McKernan, Phys. Rev. D55 2290 (1997); l.D. Lawrie, J. Phys. A21, L823 (1988); Phys. Rev. D40, 
3330 (1989); Can. J. Phys. 71, 262 (1993). 

D. Boyanovsky, I.D. Lawrie, and D.-S. Lee, Phys. Rev. D54, 4013 (1996). 
F.T. Brandt, J. Frenkel, A. Guerra, Int. J. Mod. Phys. A13, 4281 (1998). 
J. Rau and B. MuUer, Phys. Rept. 272, 1 (1996). 
X.-N. Wang, Phys. Rep. 280, 287 (1997). 

K.J. Eskola, Prog. Theor. Phys. Suppl. 129, 1 (1997); K.J. Eskola, K. Kajantie, P.V. Ruuskanen, Eur. Phys. J. CI 627 
(1998); K.J. Eskola, B. MUer, X.-N. Wang, Phys. Let t. B374, 20 (199 6 ); K. J. Eskola, N ucl. Phys. A590, 383c (1995). 



G. Aarts and J 

A. Krasnitz and R. Venugopalan |hep-ph/9909203| ; |hep-ph/9905319 ; 



Smit, Nucl. Phys. B555, 355 (199 9 ); |hep-ph/9906538|; |hep-ph/9909040 



i |hep-ph/9909040i 
" ^iep-ph/9809433t ^ 



R. Venugopalan, |hep-ph/9907209|. 

A. Rajantie and M. Hindma rsh, hep-ph/9904270 

W. Poeschl and B. Mueller, |nucl-th/9812066| ; |hcp-ph/9811441 . 

S. Jeon and L. Yaffe, Phys. R ev. D53, 5799 (1996); S. Jeon, Phys. Rev. D52, 3591 (1995). 
D. Bodeker, ^iep-ph/9903478| ; Phys. Lett. B426, 351 (1998). 
P. Arnold, D.T. Son and L.G. Yaff e, Phys. Rev. D59, 105020 (1999). 



ep-ph/9808332 tiep-ph/970632e 



J.-P. Blaizot and E. l ancu, |hcp-ph/990338£| ; |hep-ph/9906485 
M.A. Valle Basagoiti, |hep-ph/9903462|. 
D.F. Litim and C. Manuel, [hep-ph/99062l"o 
P. Danielewicz, Ann. Phys. (N.Y.) 152, 239 (1984). 

D. Boyanovsky, H.J. de Vega, R. Holman, and M. Simionato, 
Vega, Phys. Rev. D59, 105019 (1999). 

J.-P. Blaizot and E. lancu, Phys. Rev. Lett. 76, 3080 (1996); Phys. Rev. D55, 973 (1997); ibid 
Takashiba, Int. J. Mod. Phys A 11, 2309 (1996). 

R.D. Pisarski, Phys. Rev. Lett. 63, 1129 (1989); Phys. Rev. D47,5589 (1993). 

E. Braaten and R.D. Pisarski, Nucl. Phys. B337, 569 (1990); ibid. B339, 310 (1990). 

R.D. Pisarski, Physica A 158, 146 (1989); Phys. Rev. Lett. 63, 1129 (1989); Nucl. Phys. A525, 175 (1991) 
R.D. Pisarski, Nucl. Phys. B309, 476 (1988). 

N.P. Landsmann and C.G. van Weert, Phys. Rep. 145, 141 (1987). 

T. Altherr and D. Seibert, Phys. Lett. B333, 149 (1994); T. Altherr, Phys. Lett. B341, 325 (1995). 



Phys. Rev. D60, 065003 (1999); D. Boyanovsky and H.J. de 

56, 7877 (1997). K. 



Rev. Lett. 82 2840 (1999); S. Chiku and T. Hatsuda, Phys. Rev. D58 
286 (1995). 



T. Hatsuda, T. Kunihiro, and H. Shimizu, Phys. 
076001 (1998); ibid. Phys. Rev. D57, 6 (1998). 
U. Kraemmer, A. Rebhan, and H. Schulz, Ann. Phys. (N.Y.) 238, 
M.H. Thoma and C.T. Traxler, Phys. Lett. B378, 233 (1996). 

D. Boyanovsky, H.J. de Vega, R. Holman, S.P. Kumar, and R.D. Pisarski, Phys. Rev. D58, 125009 (1998). 
S. Leupold and M.H. Thoma, tiep-ph/9908460 . 
J. Schwinger, J. Math. Phys. 2, 407 (1961). 

K.T. Mahanthappa, Phys. Rev. 126, 329 (1962); P.M. Bakshi and K.T. Mahanthappa, J. Math. Phys. 41, 12 (1963). 
L.V. Keldysh, JETP 20, 1018 (1965). 

K.-C. Chou, Z.-B. Su, B.-L. Hao, and L. Yu, Phys. Rep. 118, 1 (1985). 

D. Boyanovsky, H.J. de Vega and R. Holman, in Proceedings of the Second Paris Cosmology Colloquium, Observatoire 
de Paris, 1994, edited by H.J. de Vega and N. Sanchez (World Scientific, Singapore 1995), pp. 127-215; in Advances in 
Astro fundamental Physics, Erice Chalonge School, edited by N. Sanchez and A. Zichichi (World Scientific, Singapore 1995); 
D. Boyanovsky, H.J. de Vega, R. Holman, D.-S. Lee and A. Singh, Phys. Rev. D51, 4419 (1995); D. Boyanovsky, H.J. de 
Vega, R. Holman and J. Salgado, Phys. Rev. D54, 7570 (1996). 



39 



[42] D. Boyanovsky, H.J. de Vega and R. Holman, in Current Topics in Astrofundamental Physics, Vth Erice Chalonge School, 
edited by N. Sanchez and A. Zichichi (World Scientific, Singapore 1996), pp. 183-270; D. Boyanovsky, H.J. de Vega, R. 
Holman and D.-S. Lee, Phys. Rev. D52, 6805 (1995); D. Boyanovsky, M. D'Attanasio, H.J. de Vega and R. Holman, Phys. 
Rev. D54, 1748 (1996). 

M. Le Bellac, Thermal Field Theory (Cambridge University Press, Cambridge, England, 1996). 
R.R. Parwani, Phys. Rev. D45, 4695 (1992). 
L. Dolan and R. Jackiw, Phys. Rev. D9, 332 (1974). 
P. Elmfors, K. Enqvist, and I. Vilja, Nucl. Phys. B422[FS], 521 (1994). 

L.-Y. Chen, N. Goldenfeld and Y. Oono, Phys. Rev. Lett. 73, 1311 (1994); Phys. Rev.E54, 376 (1996)^ 

T. Kunihiro, Prog. Theor. Phys. 95, 503 (1995); ibid. 97, 179 (1997); S. Ei, K. Fujii and T. Kunihiro, |hep-th/9905"088 
H.J. de Vega and J.F.J. Salgado, Phys. R ev. D56 6524, ( 1997). 



LL. Egusquiza and M.A. Valle-Basagoiti, |hep-th/9611143 
V.E. Zakharov, V.S. Lvov, G. Falkovich 



T. Kunihiro, Prog. Theor. Phys. Suppl. 131, 459 (1998). 
Kolmogorov Spectra of Turbulence I (Springer- Verlag, Berlin Heidelberg, 1992). 



Rev. 
Lee, 



D54 
Y.J. 



4066 (1996). 
Ng, and S.-Y. 



Wang, Phys. Rev. D59, 105001 (1999); S.-Y. Wang, D. 



A. Bochkarev and J. Kapusta, Phys. 
D. Boyanovsky, H.J. de Vega, D.-S. 
Boyanovsky, H.J. de Vega, D.-S. Lee, and Y.J. Ng, |hep-ph/990221E . 
R.D. Pisarski and M. Tytgat, Phys. Rev. D54 2989 (1996). 

A. Jakovac, A. Patkos, P. Petreczky, Zs. Szep, hep-ph/9905439. The authors thank A. Patkos for illuminating correspon- 
dence on these issues. 

B. D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 
D.H. Rischke, Phys. Rev. C58, 2331 (1998). 

L.P. Csernai, P.J. Ellis, S. Jeon, and J.I. Kapusta, |nucl-th/990802C . 

D. Boyanovsky, H.J. de Vega, D.-S. Lee, and S.-Y. Wang (work in progress). 

T. Altherr, Phys. Lett. B341, 325 (1995). 

L Dadic, Phys. Rev. D59, 125012 (1999). 

P.P. Bedaque, Phys. Lett. B344, 23 (1995). 

A. Niegawa, Phys. Lett. B41 6, 137 (1998). 

A. Niegawa, |hep-th/9905090 . 

C. Greiner and S. Leupold, Eur. Phys. J. C8, 517 (1999). 

M.E. Carrington, H. Defu and M. Thoma, Eur. Phys. J. C7, 347 (1999). 

D. Boyanovsky, D. Brahm, R. Holman and D.-S. Lee, Phys. Rev. D54, 1763 (1996). 

L. Faddeev and A. Slavnov, Gauge Fields: An Introduction to Quantum Theory (Addison- Wesley, Redwood City, CA, 
1991). 

K. James and P.V. Landshoff, Phys. Lett. B251, 167 (1990). 

The original method is due to Dirac: P.A.M. Dirac, Lectures in Quantum Mechanics (Yeshiva University, New York, 
1964). See also: B. Hatfield, Quantum Field Theory of Point Particles and Strings (Addison Wesley, Redwood City, 1992); 
A. Hansson, T. Regge, and C. Teitelboim, Constrained Hamiltonian Systems (Academia Nazionale del Lincei, Rome, 
1976); M. Henneaux and C. Teitelboim, Quantization of Gauge Systems (Princeton University Press, 1992), and references 
therein. 



40 




t + 




'+ k 



t + 



t' - 



(b) 






(c) 



FIG. 1. The Feynman diagrams contribute to the quantum kinetic equation for a self-interacting scalar theory up to two 
loops order. The tadpole contributions, (a) and (b), are canceled by a proper choice of A. 




(a) 



(b) 

FIG. 2. (a) The Feynman diagrams that contribute to the quantum kinetic equation for the pion distribution function. The 
solid line is the sigma meson propagator and the dashed line is the pion propagator, (b) The only contribution on-shell is the 
decay of a sigma meson into two pions minus the reverse process. 
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(b) 

FIG. 3. (a) The Feynman diagrams that contribute to the quantum kinetic equation for the sigma meson distribution 
function. The sohd hne is the sigma meson propagator and the dashed hne is the pion propagator, (b) The only contribution 
on-sheU is recombination of two pions into a sigma meson minus the decay of a sigma meson into two pions. 
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FIG. 4. The Feynman diagrams that contribute to the quantum kinetic equation for the charged scalar distribution function 
to lowest order in . The dashed and wavy lines are the HTL resummed longitudinal transverse photon propagators respectively, 
and the solid line is the HTL resummed scalar propagator, (a) Contribution from longitudinal photons, (b) Contribution from 
transverse photons. 
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